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Abstract 

We consider a quasi-linear parabolic (possibly, degenerate) equation with nonlinear dynamic boundary con- 
ditions. The corresponding class of initial and boundary value problems has already been studied previously, 
proving well-posedness of weak solutions and the existence of the global attractor, assuming that the non- 
linearities are subcritical to a given exponent. The goal of this article is to show that the previous analysis 
can be redone for supercritical nonlinearities by proving an additional Z;°°-estimate on the solutions. In 
particular, we derive new conditions which reflect an exact balance between the internal and the boundary 
mechanisms involved, even when both the nonlinear sources contribute in opposite directions. Then, we show 
how to construct a trajectory attractor for the weak solutions of the associated parabolic system, and prove 
that any solution belonging to the attractor is bounded, which implies uniqueness. Finally, we also prove for 
the (semilinear) reaction-diffusion equation with nonlinear dynamic boundary conditions, that the fractal 
dimension of the global attractor is of the order u^^-^^^^ as diffusion U — )■ 0"^, in any space dimension 



N ^ 2, improving some recent results in 23 



1. Introduction 

Let us consider the following partial differential equation 

dtu - div (a (|Vm|^) Vm) +/(«) = h{x), (1.1) 

in f2 X (0, +oo), where f2 is a bounded domain in M^, > 1, with smooth boundary F := dQ, 
a is a given nonnegative function, and / and hi are suitable functions. The mathematical 
literature regarding equation (11.11) subject to all kinds of homogeneous boundary conditions 
is fairly vast. We recall that global well-posedness results for (11.11) with Dirichlet or Neumann 

type of boundary conditions can be found in 0, 0, a a 

(see also 



addition, the analysis of dissipative dynamical systems generated by equations like (11.11) 
was carried out in a number of papers mainly devoted to the asymptotic behavior of strong 
solutions as 31,|4ll , and to establish the existence of global and/or exponential attractors 
(see, for instance, jli a, 15, 37, 46, [52[). For other classical results concerning the long term 



dynamics of ( 11. ip we also refer the reader to [a, llSl 
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All the mentioned results are mainly concerned with standard boundary conditions (that 
is, Dirichlet's and/or Neumann's). Let us now consider dynamic boundary conditions of the 
form 



dfU + b{x) a (I Vm|^) dnU + g 



u 



h2{x), 



'1.21 



on r X (0, +00), where g and /i2 are suitable functions defined on F, and b G L°° (F) , b > 
bo > 0. This type of boundary conditions arises for many known equations of mathematical 
physics. They are motivated by problems in diffusion phenomena 12 



ilJilITi, M, H H So, 



5l[ | , reaction-diffusion systems in phase-transition phenomena 
hydrodynamics 3^, l39||, models in climatology 40|, and many others. For possible physical 



26l . l47l | , special flows in 



interpretations of (11. 2p for problem (11.11) . we refer the reader to [2^] (cf. [30] also). 



Problems such as (ll.ip - (ll.2l) have already been investigated in a number of papers III.II6 



171 . |27| . l36l . |49[. Constantin and Escher deal with non-degenerate boundary value problems 
with smooth nonlinearities (in particular, (a (|,^|^) > c|^|^, with c > 0) and show 

that unique (classical) maximal solutions exist in some Bessel potential spaces 



16, 17 



Such results enable the authors to investigate other qualitative properties concerning global 
existence and blow-up phenomena (see, also [lH). These results are also improved by Meyries 



36l | , still in the non-degenerate case, by assuming more general boundary conditions and by 
requiring that / (s) /s and g (s) / s are dissipative as |s| — )■ 00. A first analysis, which aimed 
at deducing only a minimal number of assumptions on the data and nonlinearities, was done 
in ji^l by assuming that / and g are subcritical polynomial nonlinearities and by allowing 
a (s) to have a polynomial degeneracy at zero. For instance, one can take 



ais 



|(p-2)/2 



for p 2. 



;i.3) 



In particular, we proved that problem (ll.ip - (ll.2l) with a{s) as in (II. 3p . subject to square- 
integrable initial data u\t=o = uq is well-posed, and then we established the existence of a 
global attractor bounded in W^''^{Q). Well-posedness for problem (II. ip . ( II. 2p for a{s) = 
|_g|(p-2)/2^ assuming monotone functions /, g was considered in [49]. The non-degenerate case 
a{s) = u > when g = 0, is discussed in detail in 23|]. The stationary case associated with 
([□])- (USD is treated in 

It is well-known that when at least one of the source terms, the bulk nonlinear term / or 
the boundary term g is present in ( ll.ip - (ll.2p . conditions can be derived on their growth rates 
which imply either the global existence of solutions or blow-up in finite time 2l|. Namely 
in the non-degenerate case, for A, G {0, ibl} with max {A, /i} = 1, / (s) := —A s and 



g (s) := — /i \s\'^'^ ^ s, solutions of 



dtu — uAu + / 



u 



hi (x) , inQ X (0, +00) 



subject to the dynamic condition 

dtU + vbd^u + g {u) = /12 (x) , on F x (0, 00) 



;i-4) 



;i.5) 



are globally well-defined, for every given (sufficiently smooth) initial data u\t=o = Uq, if 
rir2 > 1 and Ari + /ir2 > 0. Furthermore, [21j] shows that if we further restrict the growths 



2 



of ri,r2 so that ri < (A^ + 2) / (A^ — 2) and r2 < A^/ (A^ — 2), then the global solutions are 
also bounded. On the other hand, if A = 0, /i = 1, then some solutions blowup in finite 
time with blowup occurring in the L°°-norm at a rate (t — T^,)~^^^~^^ , for some additional 
conditions on uq and r2- In the same way, when /i = and A = 1, then some solutions blowup 
in finite time with a blowup rate which depends on ri and uq (see Q). The occurrence of 
blow up phenomena is closely related to the blowup problem for the ordinary differential 
equation 

ut + h{u) = 0, (1.6) 

where either h = f or h = g. More precisely, it is easy to see that solutions of the ODE (11. 6p 
are spatially homogeneous solutions of either equation (ll.4p or equationf ll.Sp . and so if these 
solutions blowup in finite time so do the solutions of (II. 4p . (II. 5p (see (43| for further details, 
and additional references). Similar results showing the same behavior are also derived for 



the parabolic system (ll.ip - (ll.2p in 48|, stating sufficient conditions for the functions a, f 
and g so that blowup in finite time occurs in the L°°-norm. In particular, it was shown, for 
odd functions /, g and initial data Uq G L°° [Q) fliy^'^ (Q) and some additional conditions on 
Mo, that there are solutions that blow up in finite time with an upper bound on the blowup 
time which can be determined precisely. 

The main goal of this paper is to deduce more general conditions (when compared to 
conditions deduced in [ll|, |16|, ll7|, |22|, |27|, |36|, |49|, [53, on the reactive and radiation terms 



/ and g, respectively, which imply that problem (ll.ip - (ll.2p is dissipative in a suitable sense, 
and that it possesses a (possibly, finite dimensional) global attractor which characterizes 
the l ong -term behavior of the parabolic system under consideration. Recently in 4^ (see 



also 55| for some extensions), the authors have considered the semilinear parabolic equation 



(II. 4p subject to nonlinear Robin boundary conditions 

udnU + g{u) =0 on nx {0,+oo), (1.7) 

and they derived sufficient conditions on / and g, which imply dissipativity for such prob- 
lems. In particular, they have obtained a general balance between / and g, allowing for a 
real competition between both the two nonlinear mechanisms which may work in opposite 
directions, one fighting for blow-up in finite time, the other for dissipativity. Then, they also 
proved the existence of a compact attractor in (Q) , assuming that the growth of / and g 
is subcritical. Their method relies essentially on the fact that problem (II. 4p . (I1.7p possesses 
a Lyapunov functional, which can then be used to show either dissipativity of (II. 4p . (II. 7p . 
by exploiting some Poincare type inequality (see (12. 9 p below), or blow up of some solutions. 

Our goal is to extend these results in several directions, by working instead with a class 
of degenerate parabolic equations, such as (II. ip . and then by subjecting (II. ip to dynamic 
boundary conditions of the form (II. 2p . Moreover, we also wish to consider nonlinearities 
with arbitrary polynomial growth at infinity. We aim to construct weak (energy) solutions 
with the help from a different (than in i^l) approximation scheme, which is based on the 



existence of classical (smooth) solutions for a (strictly) non-degenerate system associated 
with (ll.ip - (ll.2p . Let F and G be the primitives of / and g, respectively, such that F (0) = 0, 
G (0) =0. Even though a natural energy functional exists for suitable approximates of the 
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problem f ll.ip -( !L2|) . at the moment it is not clear how to prove that this energy (see 27, 
(1.5)]), namely. 



(u) := J [a{\Vuf)\Vuf + F{u)-hi (x) u] dx + J [G {u) - h2 (x) m] 



dS 



is in fact a Lyapunov function for f ll.ip -( ITT2l) when a (s) = \s\^^ p 7^ 2, due to a lack of 
regularity of the weak solutions (see, however [23], and its references, when a{s) = u > 0). 



Therefore, the method in [43| which relies on the use of a Lyapunov function does not seem 
applicable to our situation here. Indeed, when one is dealing with gradient systems with a 
set of equilibria which is bounded in the phase space where Sfi^r (^0) < 00 , one could avoid 
to prove the existence of a bounded absorbing set and directly show the existence of the 
global attractor for subcritical nonlinearities. However, since we wish to construct global 
attractors for (ll.ip - fll.2p . under no essential growth assumptions on the nonlinearities, we 
prefer to prove the existence of a uniform dissipative estimate which can be also easily 
adapted to nonautonomous generalizations. Another difficult step that we need to overcome 
is the uniqueness problem for the weak (energy) solutions considered here. Indeed, we wish 
to deduce sufficiently general conditions on / and g without excluding the scenario based 
on which these functions are competing at infinity. When / and g are both dissipative, 
i.e., when / and g are both polynomial potentials of odd degree with a positive leading 
coefficient, uniqueness holds, and thus, the dynamical system associated with the parabolic 
system fll.ip - fll.2p can be defined in the classical sense. However, if at least one of the 
polynomial nonlinearities possesses a negative leading coefficient (for instance, suppose that 

g [s) ^ Cg\s\^'^~^ s, as |s| — )■ 00, (1.8) 

for Cg < 0), uniqueness is not know, and in this case the classical semigroup can be defined 
as a semigroup of multi-valued maps only. In order to investigate the long-term behavior 
of the degenerate parabolic system fll.ip - fll.2p . we will employ the trajectory dynamical 
approach, which allows us to avoid the use of unfriendly multivalued maps, and to apply 



the usual theory of global attractors (see, e.g., [ij] for the general theory). We strongly 
emphasize that non-uniqueness of the weak solutions constructed here is only a feature of 
the nonlinear interplay between the two nonlinear mechanisms, and is not related to the 
smoothness of the functions involved, as this is usually the case for other PDE's. See, 



e.g., [IJ, |38|, |45| and references therein, for applications for which the uniqueness is not 
yet solved, such as, hyperbolic equations with supercritical nonlinearities, reaction-diffusion 
systems, and so on. In fact, even when our functions / (s) and g (s) are assumed to be 
(locally) Lipschitz, uniqueness of the weak energy solutions is not known unless / (s) and 
g (s) are monotone increasing for all |s| > Sq (for some Sq > 0), which is quite restrictive. 
Finally, to make matters worse, for boundary nonlinearities that satisfy (11. 8p . there may be 
solutions of f ll.ip -( lL2|) which blowup in finite time at some points in Q, unless the internal 
mechanism governed by nonlinear fiux and reaction is sufficiently strong to overcome the 
boundary reaction. Therefore, it is also essential to deduce some kind of optimal general 
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conditions between the bulk and boundary nonlinearities that will only produce bounded 
(non-singular) solutions for all positive times and arbitrary initial data. The main difficulty 
here is, of course, to establish the asymptotic compactness for the system fll.ip - fll.2l) under 
some general conditions involving /, even when at least one has a bad sign at infinity, and 
to verify that any solution belonging to the attractor, is bounded. Consequently, we obtain 
the uniqueness on the trajectory attractor even for competing nonlinear mechanisms. 

To better understand the larger scope of our results, we shall illustrate their application 
to the react ion- diffusion equation (11. 4p . (II. 5p for supercritical nonlinearities, that is, for 
functions / and g satisfying the growth assumptions 

lim = (ri - 1) c/, hm = {r, - 1) c„ (1.9) 

for some arbitrary r2, ri > 1, with max (ri, > 2, and some c/, Cg G M\ {0}. Of course, our 
results below hold under more general assumptions on f,g, see Section 3. In (II. 9p . we say 
that / is dissipative if c/ > and non-dissipative if c/ < (the same applies to g). Let us 
assume bounded hi (x) and h2 (x). When both nonlinear terms cooperate, i.e., both / and 
g are dissipative and 

f (y) > -Cf, g' (y) > -c„ for all y g M (1.10) 

(for some Cf,Cg > 0), then problem (11.40 . (11.50 is well-posed and possesses a global attractor 
Agi in the classical sense, bounded in Z := W^'"^ (Q) nL°° (fi) , regardless of the size of ri and 
r2 (see Section 3.1; cf. also (23j, when g = 0). When the conditions (ll.lOp do not hold, we 
recall that uniqueness of weak solutions is not known in general. However, if / and g are still 
dissipative, we can prove that the reaction-diffusion system (II. 4p . (11.50 possesses a (strong) 
trajectory attractor Atr, which is bounded in L°° {M.^; Z) . Moreover, uniqueness holds on 
the attractor Atr so that the long-term behavior of (11.40 . (11.50 can be also characterized by 
a regular global attractor 

Agi := Atr (0) , (1.11) 

which can be defined in the usual sense of dynamical systems (cf. Section 3.2). 

For the case of competing nonlinearities, the following scenarios are possible: 
Case (i): For the case of bulk dissipation (i.e., c/ > 0) and anti-dissipative behavior at 
the boundary F (i.e., Cg < 0), problem (11.40 . (II. 5p has at least one globally-defined weak 
solution, which is bounded, if 

max (r2, 2 (r2 — 1)) < ri. 

Equality can be also allowed if the boundary condition is homogeneous, i.e., if /i2 = 0. 
Moreover, there exists a (strong) trajectory attractor Atr, bounded in L°° (IR+; Z), such that 
solutions are unique on the attractor. Thus, (ll.4p - (ll.5p also possesses the (smooth) global 
attractor Agi, defined as in (11.110 . 

Case (ii): On the other hand, in the case of boundary dissipation (cg > 0) and internal 
non-dissipation (c/ < 0), for every L^-data we obtain that, if r2 = ri = 2 (which imply that 
/ and g are sublinear) and /i2 = 0, and 

(cf + a)-' Cg)u > 2{CnCgr2y, (1.12) 
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where 



a 



{b{x))-^dS, 



and Cfi is a proper Sobolev-Poincare constant (see Section 2, assumption (H4)), then fll.4p - 
fll.Sp is well-posed in the classical sense and has a global attractor Agi, bounded in Z (in 
that case, |/'(2/)| and [(/'(i/)! are bounded for all ?/ G M by f ll.9p . see Proposition 12. 8p . 
We note that the nonlinear balance condition established in Section 3 can only be used to 
show the existence of a (strong) trajectory attractor, whenever the interior non-dissipative 
term / is sublinear as long as the boundary mechanism stays dissipative (i.e., Cg > 0), 
such that g suitably dominates / (see, ( I1.12p ). It would be interesting to see whether one 
can still construct attractors for the case of a superlinear non-dissipative function /, and 
some dissipative boundary function g of arbitrary growth. However, as we show at the end 
of Section 3, we will see that, in this case, the superlinear growth of the nondissipative 
function / produces blowup in L°°-norm of some solutions even for arbitrary nonlinearities g 
(see Section 3.4). Hence, the nonlinear boundary conditions (11. 5p for equation (11. 4p cannot 
prevent blowup of some solutions of (ll.4p -( !L5|) as long as the internal nonlinearity is strongly 
non-dissipative (for instance, when / satisfies (II. 9p with ri > 2, such that c/ < 0). Our 
main result in Section 3.4 (see. Theorem I3.28P extends some results in [3] for the parabolic 
equation (II. 4p subject to linear dynamic boundary conditions (i.e., when g = in (II. 5p ). and 



the special cases treated in [21|, |35j. In this sense, the nonlinear balance conditions derived 
in Section 3, which imply that the system (ll.4p - (ll.5p is dissipative, are optimal. 

Finally, exploiting known parabolic regularity theory for PDE's of the form (II. 40 - 01. 5p . 
the regularity of the solution for (11.40 - 01.50 increases as the functions /, g and the domain Q 



become more regular (see Remark 13.221 cf. also 23j and references therein). In particular, 
the global attractor Agi consists of (smooth) classical solutions which are defined for all 
times. Thus, using this additional regularity that solutions of (ll.4p - (ll.5p enjoy on Agi (for 
all the above cases), we obtain an explicit upper bound on the fractal dimension of Agi for 
this react ion- diffusion system by imposing weaker assumptions on the nonlinearities than in 



23(1 . In particular, for any N > 2 and for as long as (II. 9p with Cf > and Cg > holds, we 
have 

Coz/-(^-i) < dim^ Agi < Ci (1 + z/-(^-i)) , (1.13) 

for some positive constants Co,Ci which can be computed explicitly (see Section 3.3). The 
lower bound in (I1.13P was established in j23[, assuming dissipative / and (homogeneous) 
linear boundary equations (i.e., g = 0). We note that, for each fixed z/ > 0, there is a 
discrepancy between the upper and lower bounds in (11.130 as Co depends only on Q, P, /, 
g, whereas Ci is also a function of the L°°-norms of the sources hi,h2. However, we observe 
that both the upper and lower bounds are of the order u~^^~^^ as z/ — )■ 0"*", cf. (I1.13p . 
When = 1, the dimension of Agi is of the order z/~^/^, as z/ — )■ 0"*". We recall that, for 
the reaction-diffusion equation (II. 4p with the usual Dirichlet or Neumann-Robin boundary 



conditio n (|1.7D . we have upper and lower bounds of the order u for any > 1 (see. 



e.g., [42|, |45|). Thus, we have a much larger estimate (as a function of diffusion, as z/ — t- 
for the global attractor Agi in dimension > 3 (see also [23']). 
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We outline the plan of the paper, as follows. In Section [21 we introduce some notations 
and preliminary facts, then we recall how to prove the existence and stability of L^-energy 
solutions to our system fll.ip - fll.2l) . Section [3] is devoted to the existence of a bounded 
absorbing set and, then, of the global attractor Agi for unique L^-energy solutions. Then, 
we show that weak solutions possess the L^-L°° smoothing property exploiting some kind 
of iterative argument, and consequently, deduce the existence of an absorbing set in 
In the final part of Section 121 we establish the existence of the trajectory attractor Atr for 
our problem and deduce some additional properties for Atr, especially in the non-degenerate 
case when a (s) = u > 0. A blowup result for ( ll.4p -( !T3|) is also established in the case when 
g ^ 0. Finally, in the Appendix we give some auxiliary results which are essential in the 
proofs. 



2. Well-posedness in L^-space 

We use the standard notation and facts from the dynamic theory of parabolic equations 



see, for instance, 27|). The natural space for our problem is 



U2 



Si, S2 G [1, +oo] , endowed with norm 



if -Si, S2 G [1, oo), and 



\u,{x)rdx] + i^j^\u2{x)r -^^j , (2.1) 



\U\\xoo := max{||Mi||ioo(!^), ||M2||L°=(r)} 
— ||^^i||L°°(n) + ||'"2||L°°(r)- 



We agree to denote by X* the space X''''^. Identifying each function u G W^'^{Q), with the 
vector U := it is easy to see that W^'P{ri) is a dense subspace of X'' for s G [l,C)o). 

Moreover, we have 

X^ = (n,rf/i) , s G [l,+oo], 

where the measure dfi = dx\^ © dS^/b (x)|p on Q is defined for any measurable set A C by 
fi{A) = \Ann\+ S{A n r). identifying each function 9 e C [Q] with the vector 6 = Q^), 

we have that C{Q) is a dense subspace of X* for every s G [1, oo) and a closed subspace 
of X°°. In general, any vector ^ G X"* will be of the form (g^) with 6i G L** {il,dx) and 
02 G (r, dS/b {x)) , and there need not be any connection between 9i and 62. For domains 
fl with Lipschitz boundary F, recall that we have W'^'P (fl) C L^^ (fi), with = if 
pk < N, and 1 < < 00, if = pk. Moreover the trace operator Tr£,{u) := initially 
defined for u G C'^(f2), has an extension to a bounded linear operator from W'^''p{VL) into 
L'^={dri), where := if pk < N, and 1 < < 00 if = pk. For p > Nk, we have 
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W'^'^ (Q) C C^'^ (Q) , for some I. We also recall, on account of well-known generalized 
Poincare-type inequalities (see, e.g., |33|), that 

ll^lUi.p(n) '■= W^'^WiLpmf + ^ M (2-2) 

is a norm on W^'^ (Q), which is equivalent to the usual one, for any of the following expres- 
sions for /: 

I (u) := I / \u\^ —7— ] ^ ^ (u) := i lul'^ dx 



for any I < s < p. Next, for each p > 1, we let 



\u\rj 



and endow it with the norm ||-||Yfe,p given by 

It easy to see that we can identify Y'^'P with W'^^p (Q) © W'^'^^p^p (T) under this norm. More- 
over, we emphasize that V'^'^ is not a product space and since W'''P{il) ^ W^'^^P'P (F) by 
trace theory, V'^'^ is topologically isomorphic to W^'P [VL) in the obvious way. It is also im- 
mediate that V'^'^ is compactly embedded into X^, for any p > po := 2N / {N + 2) and k > 1. 
From now on, we denote by |Hlvi/fc,p(f7) and |Hli4/fc,'j(r) the norms on W'^'P (Q) and W''''^ (T) , 
respectively. Also, (■, ■)s and (■, ■)s,r stand for the usual scalar product in and L'^ (F), 

respectively. We also agree to denote by the space V^'^, and (■, •) the duality between X 
and X*, for some generic Banach space X. 

Our first goal in this paper is to give a nonlinear balance between / and g which im- 
plies dissipativity of fll.ip - fll.2p . even when both the nonlinear terms contribute in opposite 
directions. More precisely, we wish to prove the existence of (globally well-defined) weak 
solutions, provided that the nonlinearities satisfy (possibly part of) the assumptions listed 
below: 

(HI) Let b (y) := a{\yf)y, y eR^ and assume that a e C {R,R), b E (M^, M^) satisfy 
the following conditions: 

\ai\yf)\<c^il + \yr'), V|/ € M^ 

{b{y),y)^r. = a{\yf)\y\'>u\y\P, G M^ ^ 

for some constants Ci, z/ > 0. Moreover, assume that b is monotone nondecreasing, i.e., 

{a{\yif)yi-a{\y2f)y2,yi-y2)^M>0, for all yi, ?/2 e M^. (2.4) 

(H2) f,geC^ (R, R) satisfy 

lim inf/'(?/)>0, lim inf g' (y) > 0. (2.5) 

|3/|— >+oo >+oo 
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(H3) (a) f,gEC (M, M) satisfy the growth assumptions 

\fiy)\<Cf{l + \yr-'), \g iy)\ < c, {1 + \yr^') , VyER, (2.6) 

for some positive constants c/, Cg, and some ri,r2 > 1. 
(b) f,geC (R, R) satisfy 

\f{y)\<cf {l + \yr'), 
cg\yr-c<g{y)y<cg\y\''+c, WyeR, 

for some appropriate positive constants and some ri,r2 > 1. 

(H4) Let A = {J^b-^dSy^, and suppose that g E (M,M). There exists e G (0, z//g) , 
with z/ as in (12. 3p . such that, 

f{y)y + m Xy' g{y)y- \g' {y)y + g {y)\' 
hm inf , " >0, (2.7) 

\y\-^+<x> \y\ 

for some ri > p. Here 

~ ^ r Co,b(A|l]r\ if (H3a) holds 

1 %^(A|n|)-\ if (H3b) holds, ^ ' 

and Cn,6 is the best Sobolev constant in the following Poincare's inequality: 

110 - A(0/6, l)i,rlL.(f,) < Cn,t || V0||^.(^) , s>l (2.9) 



(see, e.g., |43|, Lemma 3.1]). 

We observe that condition (H4) provides an exact balance between the two nonlinear 
mechanisms. As we shall see, this balance will depend both upon the sign and growth rate 
of / and g at infinity (cf. also, (iit). 

We have the following rigorous notion of weak solution to (Il.ip -f nr2|) . with initial condition 
u (0) = Mo, similar to j27|. 

Definition 2.1. Let p e (^, +oo) n (1, +oo), and let hi (x) G L^'i (n), /i2 (x) G L''^ (F) , 
where r'- is the dual conjugate of Vi. The pair U (t) = is said to be a weak solution if 

V {t) = u (t)|p, in the trace sense, for a.e. t G (0,T) , for any T > 0, and U fulfills 

U it) G ([0, +oo) ; X^) n Wlf^ ([0, +oo) ; (V^'^)*) 

M(t)GLL([0,+oo);iyi'Mf^)), 
V it) e ([0, +oo) ; W'-^'^^^ (F)) , 

for s = min(g, r^, r'2), q := p/ (p — 1), and 

, N N N N -l^ 
fc = max(l, , ). (2.10) 

p n p r2 
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Moreover, the following identity 

{dtU, S)^, + {a {\Vu\') Vu, Va>2 + (/ (u) , a), + {g (v) /b, a^r),^^ (2.11) 

= {hu(^)2 + (^2/&,cr|r)2_r ' 
holds for all S = (/ ) G Y'^'P, a.e. t G (0, T). Finally, we have, in the space X , 

f/(0)= =:Uo, (2.12) 

where u (0) = uq almost everywhere in Q, and v (0) = Vq almost everywhere in F. Note that 
in this setting, vq need not be the trace of Mq at the boundary. 

We can cast the weak formulation (12.111) into a proper functional equation by defining 
suitable operators. To this end, let (■, ■) denote the duality between and (V-^)*. Define 
the form 

i3p(f/i, f/2) := / a (^iVuif) Vui ■ 'Vu2dx + / \ui\P~'^UiU2dx, 
Jn Jn 

for a\\Ui= C"') eYP, i = 1, 2. Note that 

BpiUi, U2) = - I div {a (|Vmi|^) Vmi) U2dx (2.13) 
Jn 

f dS f 

+ / b{x) a (jVuif) dnUiU2-r-—- + / \ui\P~'^UiU2dx. 
Jr o[x) Jq 

It follows from Lemma HIT] (see Appendix), that for each U = E Y'^, there exists Bp{U) E 
(yP)* such that 

Bp{U,W) = {Bp{U),W), (2.14) 

for every W E Y^. Hence, this relation defines an operator Bp : Y^ ^ (V^)* , which is 
bounded. Exploiting Lemma 14.11 once again, it is easy to see that Bp is monotone and 
coercive. It also follows that Bp{YP) = (V^)* (see, e.g., [^). Thus, we end up with the 
following functional form 

dtU + BpU + TiU) = gix), (2.15) 

where G {x)={'ll[l]), and the operators Bp : D (Bp) ^ X^, : D{T) C X^ ^ X^ are given, 
formally, by 

^P^-[ bix)a{\\/u\')d^u J' ^^-^^^ 

V 9{V) J 

We aim to prove some regularity results for the weak solutions constructed in Definition 



12.11 In 271], solutions were constructed with aid from a Galerkin approximation scheme 
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by imposing additional growth restrictions on the nonhnearities /, g so that they are essen- 
tially dominated by a monotone operator associated with the p-Laplacian. However, the 
application of this scheme seems to be problematic in our context since the solutions con- 
structed with Definition 12. II are much weaker than those constructed in [23, Definition 2.3]. 
Therefore, we need to rely on another scheme which is based on the construction of classical 
(smooth) solutions to the non-degenerate analogue of fll.ll) - fll.2p . One of the advantages of 
this construction is that now every weak solution can be approximated by regular ones and 
the justification of our estimates for such solutions is immediate. Thus, for each e > 0, let 
us consider the following non-degenerate parabolic problem: 

dtu - div (a, (|Vu|^) Vu) + / (n) = hi{x), in x (0, oo) , (2.17) 
dfU + b (x) (I Vu|^) dnU + g (u) = h2{x), on F x (0, oo) , 

where (s) = a{s + e) > 0, for all s G M, subject to the initial conditions 

Me(0) =Mo., iu{0)=uo,\r. (2.18) 

Let e (0) such that 



U, (0) f/ (0) = Uo in X\ 

Then, the approximate problem f l2.17p -( l2T8l) admits a unique (smooth) classical solution 
with 

Me e ([0,t,];C°° (H)) (2.19) 



for some t^> and each e > (see [ly, ll7|, lll|). Being pedants, we cannot apply the main 



results of [l6| (cf. also [l7|) directly to equations f l2.17p -f l2TT8|) since the functions a^, /, g and 
the external forces hi, h2 are not smooth enough. Moreover, the solutions constructed this 
way may only exists locally in time for some interval [0, t^,). However, taking sequences hu G 
(H) , e (r) such that hu h in L~ (Q) , and ^2 in (r) , respectively, 

and by approximating the functions a^, f, ghj smooth ones, say, in (M, M), we may apply 
Remark [3l3] below for the solutions of the approximate equations, and deduce the existence 
of a globally defined in X°°-norm solution to (12. 171) - (12. 181) . Indeed, taking advantage of the 
fact that G C°° (fi) , the global X°°-a priori bound for guarantees its global existence 
in at least V^nX°°-norm, which turns out to be sufficient for our purpose. As we shall see in 
the next section, this bound can be naturally obtained under the above assumptions on the 
nonhnearities by performing a modified Alikakos-Moser iteration argument (see Theorem 13.21 
and Remark 13.31 below). 

We shall now deduce the first result concerning the solvability of problem f ll.ip -( |L2l) . 

Theorem 2.2. Let a, f and g satisfy either the assumptions (HI), (H3a), (H4) with 

max (r2, g (r2 - 1)) < ri, (2.20) 

or (HI), (H3b) and (H4)- Then, for any initial data Uq G X^, there exists at least one 
(globally defined) weak solution U (t) in the sense of Definition \2. li 
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Proof. We divide the proof into several steps according to the different hypotheses being 
used. 

Step 1. (i) We shall now derive some basic apriori estimates for U = (t) when p < N, 
assuming that (HI), (H3a) and (H4) are satisfied. The case p > N can be treated analogously. 
The following estimates will be deduced by a formal argument, which can be justified by 
means of the approximation procedure devised above. Also, for practical purposes c will 
denote a positive constant that is independent of time, e > and initial data, but which 
only depends on the other structural parameters. Such a constant may vary even from line 
to line. Note that the smooth solutions constructed in fl2.17p - fl2.18p also satisfy the weak 
formulation (12. lip . Thus, the key choices a = Ue (t), a\r = (t) /b in (12. lip are justified. 
After standard transformations, in view of assumption (HI), we obtain 



~ + HI Vn (t)r,,(^) + (/ {u (t)) , u (t)), + {g {v (t)) , V (t) /b) 

< {hi,u{t))^ + {h2/b,v{t))^^. 



2,r 



(2.21) 



Following [43|, we can now write 

(/ (m),m)2 + {g {v),v/b)^j. 

= if (u) u + cbQ {u) u, 1)2 -CB{g{u)u- \{g (v) v/b, 1)2 p 

where 



(2.22) 



cb 



\n\ 



A 



b-^dS 



-1 



Applying inequality (12. 9p to the last term on the right-hand side of fl2.22p yields 



cb 



g{u)u-\{g (f)f/6, 1) 



2,r 



(2.23) 



<C^\\V{g {u)u) 



{epY^^ q 



Cc 



g {u)u + g{u)) Vu 



Li(n) 



g {u)u + g{u) 



Li{n) 



with Cq = Cn,fe (A ^ , and we recall that q is conjugate to p. Since assumption (H4) 
holds for some e G (0, z//g) , then from f l2.22p -f l2.23p . we obtain 

1 d 



< {hi,u{t))^ + {h2/b,v{t))^-^ + c, 



+ c \\u 



ri 



(2.24) 



for some positive constant c > 0, independent of U, t and e. Exploiting the estimate in 
Lemma [4.21 (see Appendix), and then using Holder and Young inequalities, we can bound 
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the term on the right-hand side of (12.241) by 



+ 



6 II Vnr^,(^) + Cs ( ||^||I,(^) + \\h,/b\\\ + 1 

L 2(r) 



with 7 = max(r2,p(r2 — 1) / (p — 1)), for a sufficiently small 6 > 0, and sufficiently large 
Cs > 0. Since by assumption 7 < ri, we can control the L'^'-norm of u in terms of the 
L'^^-norm of the solution (i.e., \\u\\J^j < k ||m||2Vi + C^, k ^ !)• Thus, we get for suitable 
choices of e G (0, u/q) and 6 G (0, z//g) , the following inequality 

lit 1'^^^)"^^ + I + I IIM^)lllVi(o) (2-25) 



< c 1 + ii/iiir\ + \\h 

\ L"! "L''2(r) 

for almost all t > 0. Recalling that L^^ {Q) C {Q) (we have, ri > p), we can now integrate 
this inequality over [0,T] to deduce 

U, G ([0, T] ; X^) n LP ([0, T] ; V^) , (2.26) 
u, G ([0,T] X fi), 

uniformly with respect to e > 0. On account of these bounds, we get 
Bp,, m G L'^ ([0, T] ; (V^)*) C ([0, T] ; (V'^'^') 



uniformly in e > 0, for any k > 1 (cf. Lemma [4.11 see also [27[). Here Bp,, is the monotone 
operator associated with the function a, (see (12.141) . (12.161) ). Moreover, by Lemma \^?2\ in the 
Appendix, we get at once 

V, = u,\r e L'-^ ([0, T] X F) n ([0, T] x F) , (2.27) 

uniformly in e. Due to assumption (H3a), from (I2.26l) - (l2.27p . we deduce 

G U''{[0,T] xn)x U'^{[0,T] X F). (2.28) 

Thus, J-' (Ue) is uniformly (in e) bounded in L*([0,T] ; X^), which implies 



g (x) - Bp,,U, -J'm e U ([0,T] ; (V'^'P)' 



with s = min(g, t^, rg) > 1. Therefore, dtUe is bounded in L'^ ([0, T] ; (V'^'P)*) , uniformly with 
respect to e > 0, for some k>l. Indeed, having chosen k so that W^'^ C (fi) C (fi) 
and ly'^^^/P'P c (F) (in particular, it holds N^'^ C X^^'''^ with continuous inclusion). 
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so if V" G ti^gn y ^ (yk'Py, Thus, equation (EH]) holds as an equahty in 

U ([0, T] ; iy^'P)*) and it can be considered distributionally in the space V' ([0, T] ; (V^'P)*) . 
The existence of a weak solution is then based on monotone operator arguments, followed by 
a passage to limit as e — )■ 0, and can be carried over exactly as in the proof of |27l . Theorem 
2.6]. We shall briefly describe the details below in Step 2. 

Step 1. (ii) We will now deduce the apriori bounds (I2.26p -( !2l27l) . if one assumes (H3b) 
instead of (H3a). According to (I2.2ip . in light of inequality (H3b) for g, we have the following: 

^ If/ (t) + u \\Vu (t)||i,(^) + (/ {u (t)) , u (t)), (2.29) 



2dt 



for some positive constant c. We can write, as in f l2.22p . 

= (/ (u) u + y9 («) l)^ - Y H ^ - ^ (9 (v) v/b, l)2^r 
and argue exactly as above to get the following estimate: 

^1 \\Uml + - ^) l|V« mUn) + c{\\u mi\in) + 11^ ml%ir,i>-us)) (2-30) 
< {hi, u (t))2 + {h2/b, V (t))2,r + c 

The desired control of U = (t) in fl2.26l) - fl2.27l) can be obtained immediately from a simple 
application of Holder and Young inequalities on the terms on the right-hand side of (I2.30p . 
Thus the proof is the same as in Step l-(i). 
Step 2. It is obvious that 

dtU, + Bp,,U, + TiU,) = g (2.31) 

holds as an equality in ([0,T] ; (V'^'P)*) (th is is the same as equation (12.150 with U and 
a (■) replaced by Ue and a^{-) , respectively). From the estimates that we deduced in Step 1, 

(i)-(ii), we see that there exists a subsequence {U^} = | | (still denoted by {Ue}), such 

that as e — )■ 0, 

U,^U weakly star in ([0, T] ; X^) , 

U,^U weakly in ([0, T] ; V^) , (2.32) 
dtU, dtU weakly in L' ([0,T] ; (V'^'f)*) . 

On the basis on standard interpolation and compact embedding results for vector valued 
functions (see ^^), we also have 

U,^U strongly in ([0, T] ; X^) . (2.33) 
14 



Clearly, U eC {[0,T]; (V^'f) *) . By refining in f l233|) . u, converges to u a.e. infix (0, T) and 
We converges to v a.e. in F x (0, T) , respectively. Then, by means of known results in measure 
theory (see, e.g., [i^), the continuity of /, g, and the convergence of (12.331) imply that / (m^) 
converges weakly to f (u) in U'^{Q, x (0,T)). Moreover, g (v^) converges weakly to g (v) in 
L'^2(r X (0,T)), and thus, J^(t/,) converges weakly star to (U) in ([0,T] ; (V'^'P)*). Since 
Bp^JJ^ is bounded in ([0,T] ; (V^)*), we further see that 

Bp^JJ, E weakly star in L'^ ([0, T] ; (V^)*) , (2.34) 

and thus weakly star in ([0,T] ; (V^'P)*) , since s < q and k > 1. We are now ready to 
pass to the limit as e — in equation (I2.3ip . We have 

dtU + E + J^{U) = g (2.35) 

as an equality in L' ([0,T] ; (V^'P)*). It remains to show S = BpU, which can be proved by 
a standard monotonicity argument exactly as in 27|, Theorem 2.6]. We leave the details to 
the interested reader. ■ 

The following proposition is also immediate. 

Proposition 2.3. Let a (■) satisfy (HI). In addition, assume that f,gEC (M, M) satisfy 

cilyp -c< f{y)y <C2\y\''' +c, G M, . . 

C3\yr-c<giy)y<c,\yr + c, G M, ^^'^^^ 

for some appropriate positive constants and some ri, r2 > 1 such that max (ri, r2) > p. Then, 
for any initial data Uq G X^, there exists at least one (globally defined) weak solution U (t) 
in the sense of Definition \2. ll 

Proof. In this case, both / and g are dissipative so that we do not need to exploit the 
validity of assumption (H4). Indeed, it follows from (I2.2ip that 

Ij^ + u II Vn (t)||^,(^) + ci \\u (t)ir^(^) + C3 \\v (t)irL%(r,,-M5) (2-37) 

< {hi,u{t))2 + {h2/b,v{t))^^^ + c, 

which yields the desired control of U (t) in the corresponding spaces (I2.26p - (l2.27p with relative 



ease (see 27|], for further details). Thus, the proof is the same as in Theorem 12.21 



Remark 2.4. (i)IfU (t) is a weak solution of problem ( li. jj) -( l7T^) . in the sense of Definition 
fgJl then clearly ult) eC {[0,T]; {Y'^'P)*) . Since by duality, 

for any p G {po, oo) fl (1, oo) , k > 1, and recalling that U (t) G L°° ([0, T] ; X^) , it follows 

U{t)eC^{[0,T];X') 
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(see, e.g., Il4 - Theorem 11.1.7]). Therefore the initial value U\t=o = Uq is meaningful when 
Uq G X^. Finally, we note that in general, the assumptions (H3)-(H4) alone do not ensure 
the uniqueness of the Cauchy problem U.l\) - in~^) . 

(a) We are also allowed to have equality in l{2.20\) . i.e., 7 = max(r2,g(r2 — 1)) < ri 
if h2 = 0. Indeed, this follows once again from ^2. 24^ and Lemma which allows us to 
control surface integrals in terms of volume integrals. 

Proposition 2.5. Let the assumptions of Theorem \2.2\ be satisfied. Then any weak solu- 
tion U{t) = (;;}*]) of /fO)-/ f7:^) belongs to C([0,T];X2), such that ||t/(t)||x2 «s absolutely 
continuous on [0,T] , and 

~ \\U{t) \\l = - {B, {U (t)) , U (t)) - {U (t)) -g,U (t)) , (2.38) 

for almost all t G [0,T] . We will refer to ^2. 33) as the energy identity for the parabolic 
system n~l\ )- [r^) . 

Proof. Tliis statement follows from a generalization of a known interpolation result (see, 
e.g., 14, 34, Est). Indeed, identifying the space H* = (X^)* with H = X^, we have 

V := V^'P C X^ c (V^'P)* C [Y'^'P]* =: W, 

for any k > 1. Moreover, the following inclusion E := W'^'^^ c IL"^ C E* C W also holds 
(indeed, the dual of E is the space E* = X^i'^2, and by (IXTOj) . V^'^ C E). By virtue of 
equation (12.15^ . any distributional derivative dtU (t) from V ([0,T] ; W) can be represented 
as dtU (t) = Zi (t) + Z2 (t) , where 

(t) := -Bp {U (t)) , (t) := {U (t)) + g. 

According to (12:261) - (12:281) . Zi (t) e L'i ([0, T] ; V*) = {U ([0, T] ; V))* ,q = p/{p-l), for 

U (t) e LP ([0, T];V) + {U' ([0, T] X X U'^ ([0, T] x T)) , (2.39) 

while from (ESHj), Z2 (t) G L^i ([0,T] x fi) x L^2 ([0,T] x T) , which is precisely the dual of 
the product space in (12.391) . Thus, the claim follows, for instance, from [13, Theorem II. 1.8] 
(see also [57|, Proposition 23.23]). ■ 

We will now state some results which reflect the applicability of assumption (H4) to a 
wide range of situations. In particular, it applies to the case of competing nonlinearities / 
and g, that is, nonlinearities with arbitrary polynomial growth which satisfy (113), but when 
either one exhibits a non-dissipative behavior at infinity. Recall that p G i^^, 00) fl (1, 00). 
In the interesting case of an internal dissipation mechanism, and non-dissipative boundary 
conditions, we have the following. 

Corollary 2.6. Assume that f,g E (M, M) satisfy 

lim -f-^ = (ri - 1) c; > and lim = (rs - 1) < (2.40) 

|y|->oo \s\ \y\^co \s\ 
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with Ti > p, r2 > 1. Suppose that one of the following conditions holds: 
(i) l{2.20\) holds, i.e., max (r2, q (r2 — 1)) < ri. 
(a) h2 = 0, p < r2 < q (r2 — 1) = ri and 

CfUp~'^q>C'^cl{r2Y. (2.41) 
Then, in each case the conclusion of Theorem \2.S\ applies. 

Proof. We begin by noting that (H3a) is immediately satisfied. For sufficiently large y, we 
have 

/ (y) ~ Cf \y\''^~^ y, g (y) ~ Cg \yp~^ y, 
and f {y)y ^ Cf \y\^^ , 9 iu) y ^ Cg \y\^^ . Thus, the leading terms in (12. 7p are 

Of \yr + m A)-' Cg \yr - -%r4 (^^r \y\'^''~'^ > (2-42) 

[ep)"^ q 

for some e G (0, z//g). By assumption (i), it holds 7 = max(r2,g(r2 — 1)) < ri, so the 
coefficient of the highest order term in (12.421) is cj, which is positive. If (ii) holds, it is 
obvious that r2 < max (r2 , g (r2 — 1)) = ri, so the coefficient of the highest order term in 
(^M is 

Cf ^d(r2)% 

which is positive for some e G (0, i^/g), if (I2.4ip is satisfied. Therefore, the conditions of 
Theorem 12.21 are met, and the proof is finished. ■ 

In the case of non-dissipative polynomial behavior for /, we have the following. 

Corollary 2.7. Assume h2 = and f,gEC^ (M, M) satisfy 

lim ^ = {n-l)cf <0 and \im = (r2 - 1) Cg > (2.43) 

for some p G (]|^, 00) fl (1, 00) , and let 

(c; + m Xy' Cg) up-^q > Clcl {T2r . (2.44) 

Then, the conclusion of Theorem \2.2\ applies. 

Proof. First, it is immediate that (HSb) holds with vi = r2 = p. Obviously, in this case 
r2 = Ti = q {r2 — 1) . The coefficient of the highest order term in (I2.42p is then 

{cf + m\ Xy' Cg) - -^clir2)\ 
[ep) ' q 
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which is positive for some e G (0, v/q) , if f l2.44p is satisfied. Therefore, condition (H4) holds 
true and, thus, the assumptions of Theorem 12.21 are verified. ■ 

In the case of simultaneous internal and boundary dissipation, we can prove the following 
stability result. Note that in this proposition, the uniqueness holds in the class of all solutions 
which are constructed by means of Definition 12.11 and not only for solutions which can be 
obtained as the limit, as e — 0, of the (strictly) non-degenerate parabolic system introduced 
earlier. 

Proposition 2.8. Let the assumptions of Theorem \2.2\ be satisfied, and in addition, assume 
that (H2) holds. Then, there exists a unique weak solution to problem / (i.i)) - pT^) . which 
depends continuously on the initial data in a Lipschitz way. 

Proof. Let us consider two weak solutions Ui = (J*^ ) and U2 = (J^^ ), which according to 

Proposition 12.51 belong to the spaces in (12.39^ . and set U (t) = Ui (t) — U2 (t). Thus, U (t) 
satisfies the equation 



and 

Since 

and 



dtU + {Bp (f/i (t)) - Bp {U2 {t))) + J- (f/i (t)) - J- (f/2 (t)) = 0, 
f/(0) = t/i(0)-t/2(0). 

u (t) G v ([0, T] ; VP) n {v^ ([0, T] X ^]) X ([0, T] X r)) , 



dtU it) G L-? ([0,T] ; (VP)*) + (L^i ([0,T] x fi) x L'-^ ([0,T] x T)), 
Proposition 12.51 is indeed applicable, and we have 

\j^\\U (t) \\l + {Bp {U, {t)) - Bp {U2 {t)) , U {t)) 

= -(J-(f/i(t))-^(f/2(t)),t/(t)), 

for almost all t G [0,T] . Recalling that Bp (■) G C (V^, (V^)*) is monotone and coercive (see 
fl2.14p : cf. also Appendix), we get 

IjtWUml (2.45) 
< - (/ {m (t)) - / {U2 (t)) , u {t)) - {g {m {t)) - g {U2 {t)) , u (t)) . 

Exploiting assumption (H2) (which implies, /' {y) > —Cf and g {y) > —Cg, \/y G M, for some 
Of, Cg > 0), we obtain 

±\\U{t)\\l<2{cf + Cg) \\U{t)\\l, 

for almost all t G [0,T]. Integrating this inequality over [0,T] and applying Gronwall's 
inequahty, we deduce 

\\U {t)\\l < c^' \\U {0)\\l, , (2.46) 
which yields the desired result. ■ 
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Remark 2.9. (i) Note that assumption (H2) is only required to prove uniqueness of the weak 
solution, and is usually not required for the theory of attractors. Moreover, this assumption 
is actually too restrictive so that nonlinearities that satisfy it are not allowed to carry a 
bad sign at infinity, and thus this would automatically eliminate the scenario proposed by 
the statements of Corollaries \2.6\ and \2. 7\ Indeed, a simple observation that will be made 



in Section 3 is that actually uniqueness is necessary on the attractor only, and this can be 
obtained by deducing additional regularity estimates for the solutions. This observation is 
in particular very useful if one needs to consider entropy-related nonlinearities of the form 
f (y) = log (y), for y>0, a nd f {y) = 0, for y <0, for some I > 1. ^ 

(a) Obviously, estimate (2^4^ oi/so holds if we assume that \f {y)\ and \g {y)\ are 
bounded for all y eM.. 

As an immediate consequence of the stability result just proven above, problem (11. ip - 
f ll.2p . ( I2.12P defines a dynamical system in the classical sense. 

Corollary 2.10. Let the assumptions of Proposition l278\ be satisfied. The problem U.l\) - 
( li.i^j) . I\2.1^) defines a (nonlinear) continuous semigroup S2 (t) on the phase space X^, 

^2 (t) : X2 ^ X2, 

given by 

S2{t)Uo = U{t), (2.47) 
where U (t) is the (unique) weak solution which satisfies the energy identity Ii2.38\) . 

3. Global Attractors 

3.1. Attractors for {S2 (t) , X^) revisited 

In order to study the asymptotic behavior of (ll.ip -f fL2|) . (I2.12p . we need to derive some 
additional apriori estimates for the solutions. We shall focus our study on the case p >2 only, 
since for p E ( ]|^, 2) we need to impose slightly different assumptions on the nonlinearities, 
and so we will pursue this question elsewhere. We first aim to improve some results from 



27l | for the weak solutions constructed in Definition 12. 1^ which are unique by Proposition 
12. 8[ and to show the existence of the (classical) global attractor, bounded in X°° fl V^. We 
emphasize again that all the results below hold for any p > 2. 

The next result is a direct consequence of estimate (I2.25P of Theorem 12. 2 1 (see [271 . Section 
2, Proposition 3.3], for details). 

Proposition 3.1. Let the assumptions of either Theorem \2.2\ or Proposition \2.3\ be satisfied. 
The solution semigroup {<S2 {t)}^:^^ has a -bounded absorbing set. More precisely, 

there is a positive constant Cq, depending only on the physical parameters of the problem, 
such that for any bounded subset B C X^, there exists a positive constant = (||i?||^2) 
such that 



sup 

t>t* 



\U{t)\\^. + (^j a{\Vu (s)l') \Vu (s)l' + \u {sT^ dxds 



<Co. (3.1) 
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Our next goal is to establish the existence of a bounded absorbing set in X°°, which 
has an interest on its own. The following result extends |27l T heorem 3.7] by removing the 
additional growth conditions that were imposed on f,g in |27|, Theorem 3.7, (3.17)]. 

Theorem 3.2. Let the assumptions of Proposition [XT] be satisfied. Let hi G L°° (Q), h2 G 
L°° (T), and suppose 



lim inf > 0, 



lim inf > 0. 



(3.2) 



Then, given any initial data Uq in X^, the corresponding solution U {t) of 1^2. 
belongs to X°°, for each t > 0. Moreover, there exists a positive constant Ci, independent of 
t and the initial data, and a positive constant t+ depending on t*, such that 



sup \\U 

t>t+ 



<Ci. 



(3.3) 



Proof. All the calculations below are formal. However, they can be rigorously justified by 
means of the approximation procedure devised in Section 2 (see fl2.17p - fl2.18p ). From now 
on, c will denote a positive constant that is independent of t, e, m and initial data, which 
only depends on the other structural parameters of the problem. Such a constant may vary 
even from line to line. Moreover, we shall denote by Qr {m) a monotone nondecreasing 
function in m of order r, for some nonnegative constant r, independent of m. More precisely, 
Qr (m) ~ cm"^ as m — )• +oo. 

We begin by showing that the X'^-norm of U satisfies a local recursive relation which can 
be used to perform an iterative argument. We divide the proof of (13.31) into several steps. 
Step 1 (The basic energy estimate in X™+^). We muhiply (O) by | u, m > 1, and 
integrate over Q. We obtain 



d 



(m + 1) dt 



l'"llm+i + (/(■")' I""!™ '^u)^ + m / a (|Vm|^) I Vm|^ |Mp ^ dx (3.4) 

Jn 



-I 


a 








Similarly, we multiply f]1.2p by t> 


1 d 


V 


m+1 r 

+ b{x) 
m+i,r Jy 


(m + 1) dt 


b 


= (h2 {x) , 


\v 


m— 1 \ 




/ 2,r 



im— 1 



m—1 



-.U \V 



?T1— 1 ^'5' 

+ { g[v) 



im— 1 



b x) 



(3.5) 



2,r 



Let us first observe that, in light of assumption (12. 3p . it is easy to check 

p 



m 



{\Vu\') |V 



u\ \u\ 



\m—l 



dx > um 



P 



p + m — 1 



V \u\ 



dx. 



(3.6) 
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Adding relations fl3.4p - fl3.5p . we deduce on account of the assumptions (13.21) (indeed, it holds 

f {y)y> -cis^ - C2, giy)y> -c^s^ - C4, (3.7) 

for all 1/ G M, and some Cj > 0) and an application of basic Holder and Young inequalities, 
the following inequality 



It 



m+l I 



f)+m — 1 

V \u\ p 



dx<Qi (m) (||f/||™++\ + 1) 



(3.8) 



Here the positive constant 7 ~ m~^P~'^\ p > 2, and the function Qi (m) ~ m depends also 

on the L°°-norms of hi, /12, and of bo. 

Step 2 (The local relation). Set = p'', and define 



y, it) := ^ \u (t, Or^^'"'^ dx + l^\v {t, ■)\'^"'' ^ = \\U{t) 



(3.9) 



for all k > 0. Let t, fi be two positive constants such that t — fi/rrik > 0. Their precise values 
will be chosen later. We claim that 

yk (t) < Mk {t,^i) := c {m'^y ( sup y^^i (s) + 1)"% VA; > 1, (3.10) 

s>t—fj,/mii 

where c, a are positive constants independent of k, and '■= max {zk, h} > 1 is a bounded 
sequence for all k. 

We will now prove (I3.10p when p < N. The case p > N shall require only minor ad- 
justments (in fact, in this case we can choose any arbitrary, but fixed, Ps,qs > p in the 
embedding V^'^ (fl) C XP"'*"). For each > 0, we define 



N{p + mk-l) - {N -p){l + TUk) 
N{p + mk-l)-{N -p){l + TUk-i) ' 



Sfc := 1 - Tfc. 



We aim to estimate the term on the right-hand side of (13. 8p in terms of the X^"'"™'=-i-norm 
of U. First, Holder and Sobolev inequalities (with the equivalent norm of Sobolev spaces in 
W^'P (n) C LP' (fi), Ps = pN/ (N-p)) yield 



\u 



l+nik 



dx < 



< c 



with 



(p+mfc-l)jV \ 

u\ dx j I / 1"^ 

n J \Jn 

(p+mfc-l) P 

V \u\ p dx + 



N 



n 



\u 



dx 



dx 



(3.11) 



n 



\U 



l+rrik- 



dx 



Sk ■-- 

Pi 



Sk 



[p-l) Nrrik- 



N-p {p - 1) Nnik-i + ruk + Pi 



G(0,1) 



P + ^k — 1 

p _ 1 iV - iV - p > 0, ak := ^ G [l,p] 

1 + rrik^i 
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Applying Young's inequality on the right-hand side of (13. lip , we get 



4 .,n 



V \u\ p 



\u 



dx 



max{2/i.,ctfc} 



(3.12) 



for some positive constants Tj independent of m, and where 

Zk ■= rk/ (1 - Sk) = {pi + nik+i) I {pi + nik) > 1 



is bounded for all k. Note that we can choose T2 to be some fixed positive number since 
also depends on 7^^ ~ ""^fc"^- To treat the boundary terms in (13.81) . we define for > 1, 



(N -l){p + mk-l)-iN - p) (1 + nik) 



Xfc := 1 - yk. 



{N - 1) (p + mfc - 1) - (iV - p) (1 + nik-i) ' 
On account of Holder and Sobolev inequahties (see Section 2), we obtain 



dS 



< c 



< c 



(iV-l)(p+mfc-l) \ 

v\ ^-p dS j I / If I 

p 



V \u\ 



(p+mfc-l) 



dx + I \u 



l+mfe_ 



dS 



dx 



Vk 



(3.13) 

dS^""' 



with 



N -1 



Xk 



-Xk 



(iV-1) (p-l)m,_ 



N -p {N -l){p-l) nik-i + {p-l) rrik-i + P2 ' 
P2 := {N - 1) (p - 1) - (iV - p) > 0. 



Since Xk G (0, 1), we can apply Young's inequality on the right-hand side of (I3.13p . use the 
estimate for the L^"^"^*-' (f2)-norm of u from (I3.12p in order to deduce the following estimate: 



(3.14) 



4 Jn 



p+m^ —1 



V \u\ p 



dx + {nik) 



l+nik- 



^ dx 



max{Jfc,afc} 



for some positive constants t^, depending on ti, r2, but which are independent of m. The 
sequence 

Vk {p - 1) nik+i + PP2 



Ik '■- 



> 1 



(1-Xfe) {p~l)mk+pp2 

is bounded for all A; > 1. Inserting estimates (l3.12p -( l3TT4l) on the right-hand side of (13. 8p . 
we obtain the following inequality: 



dtyk + Ik 



V \u\ p 



dx<c{muY' iyk^i + ir 
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(3.15) 



for some positive constant ai that depends on Xj; we recall that = max {z^, 4, a^} > 1, 
and 7fc ~ rn^]^^ . 

Let now ^ (s) be a positive function ^ : IR+ — > [0, 1] such that C,{s) = for s G 
[0,t — /i/rfc], ^ (s) = 1 if s G [t,+oo) and < rrik/fi, if s G {t — fi/rk,t). We de- 

fine Zk (s) = C (s) yk (s) and notice that 

^z, is) < c (s) is) + "^y, {s) (3.16) 

as as jj, 

= C is) ^yu is) + Qi im,) (^^ rfx + ^ Ivl'^"'^ . 

The last two integrals in fl3.16p can be estimated as in f l3.12p and fl3.14p . Combining the 
above estimates and the fact that Zk < yk, we deduce the following inequality: 

-^Zk is) + crrikZk (s) < (t,/i) , for all s E [t - fx/vk, +oo) . (3.17) 

Note that c = c (/i) ~ as — )■ 0, and c (yu) is bounded if fi is bounded away from zero. 
Integrating fl3.17p with respect to s from t — fi/vk to t, and taking into account the fact that 
Zk it — f^/rk) = 0, we obtain that (t) = Zk it) < Mk it,fi) (1 — e~'^^), which proves the 
claim fIXTU]) . 

Step 3 (The iterative argument). Let now r > r > be given with r = as in (13.11) . and 
define /x = pir — r) > 1, to = r = + 1 and tk = tfc-i — fJ^/mk, k > 1. Using (I3.10p . we 
have 

sup yk it) < c ifi) (mfc)" (sup yk-i is) + 1)"\ A; > 1. (3.18) 

t>tk-l S>tk 

Note that from (13. ip . we have 

sup iyois) + l)<Co + l=:C, (3.19) 

S>tl=T 

and c = c (/i) is bounded away from zero. Thus, we can iterate in (I3.18P with respect to 
k > 1 and obtain that 

sup yk it) < icml) {cmUY" {cmUY''''-' ■ ... ■ (cm^)"'="— {C)^- (3.20) 

where ^k '■= nkUk-i-.-Uo, and 

Ak := I + Uk + UkUk-i + ... + nkUk-i-.-no, (3.21) 
Bk:= k + Uk (A; - 1) + n^nfc^i (A; - 2) + ... + nkUk-i-.-UQ. (3.22) 

Without loss of generality, let us assume that Zk > h > 1, for each k. Then, = Zk, and 
Cfc = ipi + ^k+i) I (pi + w^o)- The argument below also applies to the case when Uk = h- 
Thus, we have 

Ak < iPi + rrik) y] ■ and Bk < (pi + ruk) ■ . (3.23) 

Pi + rrii ^ pi + rrii 

1=1 1=1 
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Therefore, since 

sup yk {t) < sup yu (t) < c^'=p"^'= iff' (3.24) 

t>to t>tk-i 

and the series in (13.231) are convergent, we can take the 1 + mfc-root on both sides of (I3.24p 
and let /c ^ +CX3. We deduce 

sup \\U{t)\\^^< hm sup(3^fe(t))'/^'+"^^^ <Ci, (3.25) 

t>to=T' ^~'+°° t>t0 

for some positive constant Ci independent of t, k, U, e and initial data. The proof of Theorem 
13.21 is now complete. ■ 

Remark 3.3. (i) We can easily modify our argument in the proof of Ii3.3\) in order to show 
that the li°°-norm of the solution U (t) stays bounded for all time t > 0, if Uq is bounded in 
the X°°-norm. It suffices to note that in place of the inequality Ii3.10\) . we may use instead 
the inequality 

ykit) < Q(||f/o|lx-,supMfc(t,/i)), 

t>o 

which is an immediate consequence of Ii3.15\) . Then arguing as in the proof of Theorem \3.2l 
we also have the estimate: 

sup||t/(t)||x^ <Q(||[/o||x^,sup||f/(t)||^0, (3.26) 

t>0 t>0 

for some positive monotone nondecreasing (in each of its variables) function Q : — > M_|_ 
independent of e. 

(a) By slightly refining the arguments in the proof of Theorem \3.^ (in Step 3), it is also 
easy to show that, for each r > 0, 

sup \\U{t)\\^^<Q{r-\snv\\U{t)\\^,). 

t>2T t>T 

Arguing exactly as in the proof of Theorem 12.21 we can also obtain the following general 
balance condition between the functions f,g, implying boundedness of the solution. 

Proposition 3.4. Let the assumptions of Theorem \2.2\ be satisfied, and let hi G (Q), 
/i2 G (r). Suppose that there exist r > and yo > 0, such that for any m > 1 and 
\y\ > yo, it holds 

f (y) lyr' y + m A)"' g (y) \yr' y - ^""Z/T lyr' g' [y) v + mg (y) ' (3.27) 

[ep) ' q 



> -Qr (m) 



m+l 



for some e G (0, and some positive function : IR+ — )■ ]R_|_, (m) ~ m^, as m ^ oo. 
Then, the same conclusion of Theorem \3.S\ applies to any weak solution of problem ( fj. ijj - 
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Proof. Let us return to equations (13. 4p - (13.61) . We can write 



\v\ V 



(3.28) 



2,r 



I im— 1 \ 

[u],\u\ u). 



= (/ (u) + csg I 
- cb {g (u) u - A (^f (u) u ImT""^ , l/h)^ .^^ , 

where cb and A are as in Theorem 12.21 Applying the Poincare's inequahty (12. 9p . we have 

(g (u) u ImT"^ - a (^f (u) u , l/b)^ Jj^ 

im— 1^ 



Cb 



<C\\V{g iu)u\ur~^)\ 



C 
C 



^Vuig {u)u + mg{u) 



Ill — ± \ lit — i. / / 

u\ f Vm I \u\ 9 \g {u)u + mg{u) 



dx. 



On account of standard Holder and Young inequalities, we can estimate the term on the 
right-hand side in terms of 



1/p 



\u\ 



g {u)u + mg [u] 



dx 



(3.29) 



C 



P 



p + m — 1 



m 



p+m — 1 

V \u\ p 



< em 



im— 1 



P 

dx j 

g (u) u + mg [u) dx \ m^^^^ 



dx 



( " Yf 


p+m — 1 


V \u\ p 


\p + m-lj 





+ 



Cim-i/p 



\u\ 



im— 1 



g {u)u + mg {u] 



{epY^^ q Jn 

Recalling (I3.28p . on the basis of (13.290 . we can estimate 



dx. 



{fiu),\ur u)^ + (g{v), 



I I'm— 1 

\v\ V 



(3.30) 



> (f (u) + CBg{u) ,\u\"^ M>2 



2,r 



— em 



P 



p + m — 1 



jj-\~7n—l 

V \u\ p 



{epY^^ q 
p 

dx. 



g {u)u + mg {u) 



\u\ 



im— 1 
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Therefore, combining (13 ■4p with (13. Sp . then using (13. 6 p and (I3.27p . we arrive at the following 
inequality 



^ \\Uh^+i+i{jy -e) 



p-j-m — 1 

V \u\ p 



rfx < g,(m) (||f/||™+i + 1) , (3.31) 



for some e G (0,z//g). From this point on, the proof goes on exactly as in the proof of 
Theorem 13.21 (cf . Steps 2 and 3) . We omit the details. ■ 

We will now verify the hypothesis in Proposition 13.41 for functions that satisfy the as- 
sumptions of Corollaries 12.61 and 12.71 

Corollary 3.5. Assume that (HI) holds, and the functions f,g & (M, M) satisfy all the 
assumptions of Corollary \2.6\ -(i). Then, for any initial data Uq in X^, the corresponding 
solution U (t) of U.l\) - in~^) . li2.1S\) belongs to X°°, for each t > 0, and estimate liS. 3\) holds. 

Proof. As in the proof of Corollary 12. 6[ the leading terms on the left-hand side of (I3.27P 
are, for sufficiently large |?/| S> 1 and any m > 1, 

Cf lyr"-""-' + m X)-' c, - ^^^^^cl {r, + m- If , (3.32) 

[ep) ' q 

for some e G (0, z//g). From Corollary 12. 61 (1). it holds 7 = max{r2,q {r2 — 1)) < ri, so the 
coefficient of the highest order term in (I3.32p . for any m > 1, is c/ > 0. Therefore, the 
desired claim follows immediately from (I3.27p . ■ 

Corollary 3.6. Let /i2 = and assume f,gEC^ (M, M) satisfy 

lim / (s) = (ri — 1) c/ < and lim g (s) = (r2 — 1) Cg > 0. 

\y\-^oo |y|^oo 

Then, the conclusion of Theorem \3.S\ applies. 

Proof. The proof follows, for instance, from Theorem 13.21 since the functions /, g satisfy 
O. ■ 

Having established that the weak solution is bounded for any positive times, we also have 
the following. 

Proposition 3.7. Let the assumptions of either Theorem \3.2\ or Proposition \3.4\ be satisfied. 
Then, any solution U (t) of U.1\) - I[T^) . h2.12\) belongs to V^, for each t > 0, and the following 
estimate holds: 

sup (\\U(t)\%v + I \\dtU{s)\\l,d^ < C2, (3.33) 
t>ti \ Jt J 

for some positive constant C2, independent oft, e and initial data. 
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Proof. It suffices to show (I3.33p . We ffist recall that, using assumption (12. 3 p and the fact 
that Ttd : W^'P (^]) W^i-i/p-p (r) is a bounded map, from (El]) we get 



sup 

t>t# Jt 



t+i 



\U{s)\\l, + \\u{s^^^^^ 



L'-i(n) 



ds < a 



0' 



(3.34) 



for some positive constant Cq independent of time and initial data. 

Let us now multiply equation (II. ip by dtu {t) , then integrate over f2, and multiply equa- 
tion (II. 2p by dtv (t) /h (x) and integrate over V. Adding the relations that we obtain, we 
deduce after standard transformations, 



1 d 
2~dt 



{A{\Vu (t)n, 1), + 2{F{u (t) , 1))^ + 2{Giv (t)) , 
h2 



(3.35) 



2,r 



-2(/ii,n(t)) 



2-2(y,t;(t) 



- Wdfvl 



2,r 



lL2(n) ll^t'^llL2(r,d5/b) ' 

for all t>tj^ (with as in (13. 3p ). where we have set 

,|2 



A{\y?) 
Next, let us define 



a (s) ds, F (y) 



f{s)ds,G{y) 



g{s) ds. 



£ it) := {A{\Vu {t)f), + 2{F{u it) , 1))^ + 2{G{v (t)) , 

-2{h,,u (t))2 - 2 {h2/b, V (t))2 p + Cf,G- 



(3.36) 



Here the constant Cf,g > is taken large enough in order to ensure that S (t) is nonnegative 
(recall that F (u) and G (v) are both bounded by (13.31) ). On the other hand, on account of 
(13. 3p . one can easily check, using the fact 

(A(|Vn|'),l)2>Cp||Vn||^,(^), cp>0, 

that there exists a positive constant c, independent of t and the initial data, such that 

\\Vuit)\\%-c<£{t), (3.37) 

for t > max t*} . From (13.351) . we have 

dS{t) 



dt 



+ 2 \\dtv\\L2(j:^ds/h) = 0' ^ max{t+,t*} 



(3.38) 



Then, exploiting estimates (13.31) . (I3.34p and (I3.37p . we can apply to (I3.38P the uniform 
Gronwall's lemma (see, e.g., |45|) and find a time ti > 1, depending on t+, i:*, such that 



^ < r 
Lp{n) — 



Wt > ti, 



(3.39) 
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for some positive constant c. Summing up, we conclude by observing that f l3.33p follows from 
fl3.39p and the boundedness of the trace map Ttd : W^'^ {Vt) — H^^^^/P'P (F) . The proof is 
finished. ■ 



Finally, the above dissipative estimates and standard compactness results (see, e.g., |27|, 
Section 2]), allow us to conclude the following. 

Theorem 3.8. Let the assumptions of either Theorem \ 3.2\ or Proposition \3.4\ be satisfied 
and, in addition, let (H2) hold. Then, the dynamical system (iS2(t),X^) generated by the 
initial value problem M.1\) - [T7B^} . h2.12\} possesses the global attractor Agi C X^, which is a 
bounded subset ofY^ fl X°°. Moreover, 

lim distxn.s2 (^2 (t) B, Agi) = 0, (3.40) 

for any finite Si, S2 > 2, for all bounded subsets B o/X^. 

3.2. Trajectory dynamical systems 

In the final part of this section, we shall devote our attention to constructing the "usual" 
weak trajectory attractor and verify (using the maximum principle established in Theorem 
13. 2p that any solution, belonging to the attractor, is bounded so that uniqueness holds 
on the attractor. We will employ a slightly different construction (compared to e.g., 
and references therein) of the trajectory attractor, which also looks more natural from the 
physical point of view. Namely any weak solution U (t) of f lLip -f lTT^ . fl2.12p is included in 
the trajectory phase-space of the problem if and only if it can be obtained in the limit, as 
e — J- 0, of the corresponding solutions (t) of the approximate system fl2.17p -f l2T8|) . 

In order to define the trajectory dynamical system for weak solutions without uniqueness 
we need to introduce first the appropriate functional framework. First, let us recall estimates 
f l2.25p and f l2.30p which hold for any smooth solution U = oi the approximate problem 
f l2.17p -( pT8|) . By a standard application of Gronwall's inequality (see, e.g., (45[), we get the 
well-known estimate 

\\Uml + f Wll^. + lk(r)irL.,(^) + c|mr)||-,(r)) dr (3.41) 

< \\U{s)\\l2 e-<'-^^ + c{l+ \\h,\{\ + \\h2\\'\^^^ (1 - e-(*-^)) , 

for all t > s > 0, and some appropriate positive constant c. Here and below c = 1, if (H3b) 
is assumed and c = 0, when (H3a) holds. Let 0^'^°*^ denote the local weak topology in the 
space 

L°° (M+; X^) n (M+; V^) H (L'^^ (M+; ((])) x (M+; (F))) . 

By definition, a sequence f/„ (t) — )■ f/ (t) , as n — )■ 00, in the topology of 0^|!''°^ if, for every 
T > 0, 

Un (t) ^U{t) * -weakly in L~ ([0, T] ; X^) , 
Un (t) U (t) weakly in ([0, T] ; V^) , 

Un (t) ^ U (t) weakly in ([0, T] ; {Q)) x ([0, T] ; (F)) . 
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We recall that GXl ' is a Hausdorff and Frechet-Urysolin space with a countable topology 



space (see, e.g., j4J]). Next, let 9^ be the Banach space defined as 

:= (M+; X^) n (M+; V^) H {V^ (M+ x ^]) x L'^^ (R+ x T)) . 

Note that the unit ball of is compact in the local weak topology of 0^''°*^ (see [4 

Definition 3.9. A function U e Q\ is a solution of (EJlj-lElj, HJD with Uq G X^ if 
it solves ( fi.i)) -( l7T^) in the sense of Definition \2. il and i/iere exists a sequence e„ — )■ 0, a 
sequence Uq. G C°° (fi) , f/o. = ( ) anc? a sequence U. of classical solutions of the 

approximate problem ^2. Ill )- f2. 180 with e = e„ such that 

Uo = Xl- lim Uoe„ and U = 6+''°" - lim U,„. (3.42) 

n^oo n— >oo 

Note that fl3.42p implies in a standard way the weak convergence of 

i3p,e„f/.„ -> B,U in LL (M+; (V^)*) , (3.43) 
and consequently, the weak-star convergence of dtU^^ — dtU in 

LL (V^)*) + (Li (M+; L-i (fi)) x Li(M+; L^^ (F))). (3.44) 

This gives the strong convergence f/e„ (t) ^ f/ (t) in Cioc (V'^'P)*) (see Section 2). Thus, 
any solution U of fll.ll) - fll.2l) . (12.1 2p is weakly continuous with values in X^ (see Remark [23]), 
and for any t > 0, we have the weak-convergence 

Ue,. (t) ^ U (t) (3.45) 

in the space X^. It is important that we do not require the strong convergence in (13.45^ even 
for t = 0. 

We can now summarize the results in Section 2 by stating the following. 

Proposition 3.10. Let the assumptions of Theorem \2.2\ hold. Then, for every Uq G X^, 
there exists at least one globally defined weak solution U {t) , t G M^, of the degenerate 
problem M.l\) - n~^) . with U (0) = [/q, which can be obtained as a weak limit l[3.4^ of the 
corresponding solutions Ue„ (t) of the approximate non- degenerate parabolic system l[2.11^ - 

In order to construct the global attractor for the dynamical system associated with the 
degenerate parabolic system fll.ip - fll.2l) . we need the following definition. 

Definition 3.11. Let 

Mu (t) := inf (liminf ||f/,„ {t)\\^, : U = Ql''"' - lim f/,„, f/o = X^ - lim f/,„ (0)1 , 

where the external infimum is taken over all possible sequences of solutions of the approximate 
problem { 2.17 )- f2.18\} . which converges as e„ — t- to the given solution U of the limit problem 
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Some simple properties of this M-functional are stated below. 

Proposition 3.12. Let U be a solution of problem 111. 1\) - [T7^) . ^2.12^) and let Mjj {t) be the 

associated functional. Then, 

(a) II W 11x2 < Mu it) , for all t e M+. 

(b ) The following estimate holds: 

(Mu {t)f + 1^ (||f/ Wfv. + \\u (r)||-,(^) + c lit; (r)||-,(r)) dr (3.46) 
< {Mu (5))'e-^(*-^) +c(l + ||/ii|r\ + \\h2p, e-<'-'^) , 

where t > s > 0. 

(c) MT(h)u (t) <Mu{t + h), for all h>0, where (T (h) U) {t) := U {t + h) . 

Proof, (a) is immediate since the norm ||-||x2 is weakly lower semicontinuous and the 
convergence of U^^ to U in 0^''°^ implies the weak convergence Ue„ if) — )■ U (t) for every t. 
To prove (b), we note that due to the energy estimates for the approximate parabolic system 
driZD-dllHD (cf. Section 2), we have 



\\Ue. ml + {\\U.. Wfv. + he. (r)||lV,(^) +c\\v,„ (r)|ll%(r)) dr (3.47) 
< {s)\\l e-(*-) + c(l + \\h,rl + \\h,p, ) (1 - e-(*-)) , 

\ L i L 2{T) J 

for every U^^. By definition of M[/, for every 5 > 0, we can find an approximating sequence 
U^^ such that 

liminf ||?7,„ mW^<Mu{t)+5. 



Passing to the limit, as n — )■ cxd, in (I3.47p . we have 

[Mu m" + f {\\U{r)n. + \\u (r)||2V ,(^) + c\\v (r)||-,(,)) dr 

< liminf (^||f/,„ mil + f {\\Ue„ ir)K. + lke„ {r)\\l\(^^^ + c\\v,„ W||l%(r)) dr 

< {Mu (s) + 5f e-<'~'^ + c(l + \\hif\ + \\h2p, e-^(*-^)) , 

and since (5 — )■ is arbitrary, we get the desired inequality. The third assertion is also 
immediate since the infimum in the definition of MT(h)u (t) is taken over the larger set of 
admissible approximating sequences than the infimum in the definition of Mu(t+h)- ■ 

We are now ready to construct the trajectory-phase space, the trajectory semigroup and 
the kernel associated with the degenerate problem (ll.ll) - fll.2p . To this end, let /C+ C 9i be 
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the set of all solutions of fll.ip -( lL2|) . in the sense of Definition 12. which corresponds to all 
f/o G X2, and let 

T{h) : /C+ ^ /C+, h> 0, (3.48) 

be 

{T{h)U){t) ■.= U{t + h). (3.49) 

We shall refer to /C+ and T (h) : /C+ — )■ /C+ as the trajectory phase space and the trajectory 
dynamical system, respectively, associated with the degenerate parabolic system fll.ll) - fll.2p . 
In addition, we endow the set /C+ with the topology induced by the embedding /C+ C 9!|!''°'^ 
and we will say that a set i? C /C+ is M-bounded if 

Mb (0) := sup Mu (0) < oo. 

Note that any M-bounded set B C /C+ is bounded in the norm of G!^. Finally, a kernel 

:= (M; X^) n L'^ (M; V^) n {U' (R x Q) x (R x T)) 

consists of all complete (defined for all t G M) bounded solutions of fll.ip - fll.2l) . fl2.12p which 
can be obtained as the weak limit, as e„ — j- 0, of the appropriate solutions of the approximate 
non-degenerate parabolic system fl2.17p - fl2.18p . Namely, f/ G /C if and only if there exists 
a sequence e„ — )■ 0, a sequence of times tn — > — oo, and a bounded sequence of initial data 
Moe„ G {Vl) , ||f/oe„|lx2 — ^i ^uch that the corresponding solutions Ue„ of f l2.17p on the 
interval [t„, +oo) with initial data U^^ (t„) = Uq^^ converges weakly in 6^ to the complete 
solution U considered. 

We now recall the definition of the global attractor for the trajectory dynamical system 



{T (h) ,/C+) (see for more details; cf. also |38l. l56|) 



Definition 3.13. A set Atr C /C+ is a (weak) trajectory attractor associated with the de- 
generate parabolic system U.l\) - in~^) (= global attractor for the trajectory dynamical system 
(T {h) , /C+) ) if the following conditions are satisfied: 

(i) Atr is compact in /C+ and is M-bounded; 

(a) It is strictly invariant, i.e., T {h) Atr = Atr, h > 0; 

(Hi) It attracts the images of bounded (M-bounded) sets as h — )■ oo, i.e., for every B 
bounded in and every neighborhood 0{Atr) of Atr (in the topology ofQ^^°'^), there exists 
ho = ho {B, O) such that T{h)B C O (Atr) , V/i > /iq. 

The next theorem can be considered as the second main result of this section. 

Theorem 3.14. Let p G (jf^jOo) fl (l,oo). Let all the assumptions of Theorem \2.^ be 
satisfied. Then, the degenerate parabolic problem hl.l\) - [T7^) . ^2.12) possesses a trajectory 
attractor Atr C 6^ and the following description holds: 

Atr = nt>o (/C) . (3.50) 

Here and below, Ht^i (/) denotes the restriction on I of a function f defined on M. 
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Proof. According to general theory IJ, |38|, |56|, we are only required to check that the 



trajectory dynamical system is continuous and that it possesses a compact and M-bounded 
absorbing set. The continuity is immediate since T (h) are continuous on 0!ji''°'^. The estimate 
(b) in Proposition 13.121 guarantees that the set 

B:= \u e}C+:{Mu{0)f <2c(l + \\hi\\"\ +\\h2p, )] (3.51) 

will be absorbing for the semigroup T (h) : /C+ — > /C-|_. Moreover, this set is semi- invariant. 
This follows from Proposition 13.121 (c) since 

<{Mum\-^^'-'^ + c(l + \\h\\%+\\h,\\% ) (1 



g-c(t-^)^ 



<2c( 1+ Whip, + \\h 



,, 2|| ^' J , 

L^'i L''2(r) 

for all U E B. Therefore, T {h) B G B. It remains to show that B is compact. Due to the 
inequality (b) in Proposition I3.12[ the set B is bounded in 0*^, and therefore precompact 
in 9^''°'^. Thus, we only need to show that B is sequentially closed, i.e., if Un E B and 
U = 6^' — lim„^oo Un, then U E B as well. 

For Un G B, Mu„ (0) is bounded. By estimate (b) in Proposition [3?T2]the sequence Mu„ (t) 
is bounded for all t > 0. Moreover, since every Un is a solution of the degenerate problem 
f ll.ip -f lTT^ . there exists a sequence Un,ei of solutions to the approximate non-degenerate 
problem fl2.17p -( l2rT8|) such that ei = en,i — )■ as /c — )■ oo, and 



Without loss of generality, we may assume that 



\Mu„ (0) - ||f/n,e, mir^ ei< -, (3.52) 



for all / G N, and we may also suppose that 

liminfMc;„ (0) = lim Mu„ (0) 



n— >oo n—¥oo 



(we may pass to a subsequence in n if necessary). It remains to show that we can extract 
from {Un,ei}ni^N ^ parametric sequence which will converge to the limit function U. To 



this end, recall that the topology of 0^'^°*^ is metrizable on every bounded set of 0^ (see 44|). 
Let 6 > he such that all Un belong to the closed ball Bs of 0^. Evidently, U E Bs C B25, 
and we may also assume that Un,ei G B2S, for all n,/ G N. Indeed, the sequence Un,ei (0) 
is uniformly bounded in n, Z, due to (13.521) . and since Mjj^ (0) is bounded, then recalling 
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estimate fl3.47p . we also get that Un,ei (t) is uniformly bounded with respect to n,l and t. 
Let d (-, ■) be a metric on B2S- Thus, 

lim d {Un, U) = 0, lim d {Un, Un,,,) = 0, 

ra— i>oo fc^oo 

for every n. Therefore, for any n, there exists Iq = Iq [n) such that d [U, Un,ei) < I/'t-, for all 
/ > Iq. Thus, we have d{U, Un,ei^) — !■ as n — > oo, and therefore 

^ ra-i-oo 

Moreover, thanks to fl3.52p . e/p — )■ as — )■ oo and so [/ is a solution of the degenerate 
problem, and 

Mu (0) < liminf ||f/„,,^ (0)||^, = liminfM^/^ (0) . 

n— >oo ra— >oo 

Thus, B is indeed a compact semi-invariant absorbing set for (T (/i) , /C+) , and the desired 
attractor can now be found in a standard way, as the w-limit set of B : 



Ar = uj{B) = n T{h)B. 



h>0 

The description fl3.50p is a standard corollary of this explicit formula and the diagonalization 
procedure described above. This completes the proof. ■ 

What is the connection between the dynamical system {S2 {t) , X^) introduced at the end 
of Section 2, and the trajectory dynamical system (T (h) , /C+) constructed here? It turns out 
that, under the assumptions of Proposition I2l8| the solution U{t) of the system fll.ip - fll.2p is 
unique and, consequently, this parabolic system generates a semigroup in the classical phase 
space X^, in a standard way by the formula fl2.47p . If we define the map 

^o:/C+^X^ Uo{U{t)):=U{0), 

we see that the map IIq is one-to-one and, in fact, Hq defines a Lipschitz homeomorphism 
between /C+ and X^ (i.e., X^ endowed with the weak topology). Therefore, when uniqueness 
holds (for instance, if we require that the functions f,g satisfy (112)), 

S2{t) = UoT{t) 

the trajectory dynamical system (T (t) , /C4.) is conjugated to the classical dynamical system 
^2 (t) defined on the phase space X^ endowed with the weak topology. However, we note that 
assumption (H2) is quite restrictive as it does not allow for a competing scenario between the 
nonlinearities f,g, as proposed in e.g.. Corollaries 12.61 and 12.71 So without assumption (H2) 
the uniqueness problem for fll.ll) - fll.2p is not known, and the classical semigroup ^2 (t) can 



be defined as a semigroup of multi- valued maps only (see, e.g., [38j for further details). The 
trajectory dynamical approach allows us to avoid the use of multivalued maps and to apply 
the usual theory of global attractors to investigate the long term behavior of the degenerate 
parabohc system ([III]) -([112]), I^J^. 

As a consequence of Theorem I3.14[ we can also state the following. 
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Corollary 3.15. Let the assumptions of Theorem 3.14 hold and let B C /C_|_ be an arbitrary 
M -bounded set. Then, for every T G M+ and every s G {0,k), the following convergence 
holds: 

lim distc((h,T+h);(y^'py) {B\{h,T+h), Atr\(h,T+h)) = 0. (3.53) 

Proof. Indeed, from the fact that every U G /C+ is a weak solution in the sense of Definition 
12.11 (cf. (12. lip ), we can express and estimate (as in Section 2) the time derivative of U, i.e., 
dtU G Ll^ (M+; {Y'^'P)*) . Next, since the embedding 

[U:Ue X^) , d,U G {Y^'P)*) } C (M+; (V^'^)*) , (3.54) 

is compact for every < s < A;, then /C+ C Cioc(M+', (V^'P)*) is also compact in the sense 
that every M-bounded subset of /C+ is a precompact set in the set Cioc (V*'^)*) . Thus, 
the above convergence (I3.53P is an immediate corollary of Definition 13.131 (iii). ■ 

In the sequel, we shall also verify, under additional assumptions on the nonlinearities 
(which still allow for a competing behavior between / and g), that every solution U of the 
degenerate parabolic system is uniformly bounded on the (weak) trajectory attr actor Atr- 
Then, using this fact we can establish that the solution of (ll.ip -f lL2|) . (I2.12p is unique on the 
attractor. From now, we will always assume that p >2. 



Theorem 3.16. Let the assumptions of either Theorem \3.S\ or Proposition\3.4\ be satisfied. 



Then, for every complete solution U G JC, we have U G L°° (M; fl V^) and the following 
estimate holds: 

11^ (^)llx°°nvp — Qi^ + ll^illL°=(f7) + ll^2|lioo(r)), (3.55) 
for all t G M, for some monotone non- decreasing function Q independent of U , t and initial 
data. Moreover, U G Cioc (IR; X^) , for every U & IC. 

Proof. The proof is essentially based on the maximum principle deduced in Section 2, and 
the description of K, from Theorem I3.14[ Let f/ G /C be an arbitrary complete solution, i.e., 
let e„ — )■ 0, t„ — 7- — oo and U^^ (t) , t > tn be the sequence of solutions of the approximate 
parabohc system fl2.17p with U^^ (t„) = Uq^^, where ||f/oe„|lx2 ^ uniformly with respect to 
n. Let us now fix an arbitrary T G M. Then from the convergence Ue„ — ?■ U, we also know 
that 

f/,„ (T) ^ f/ (T) (3.56) 

strongly in X^ (passing to a subsequence in n, if necessary), since the embedding /C+ C 
Lf^^ (]R+; X^) is compact (i.e., any M-bounded subset of /C+ is precompact in Lf^^ (M+; X^)). 
It follows from Theorem 13.141 and estimate (13.471) that 



\U.„ (T)ll^. + (||f/,„ (r)ll^, + (r)||2.,(^) +c\\v,„ (r)|II%(r)) dr (3.57) 



<2c l + \\hi\\'\ +\\h2\\', 

L^'i L''2(r) 
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where the constant c > is independent of e„ and T. By modifying the proof of Theorem 
I3.2l in a suitable way, we arrive at the following inequality for the approximate solutions U^^, 

sup ||t/,„ {t)\\^^ < g(c + l|f/||^..(i,,x2)), (3.58) 
t>t„ 

< Q{1 + ||/?.i||ioo(f^) + ||^2|lLoo(r)), 

for some monotone function Q independent of e„, T, t. Passing to the limit, in a standard 
way in (13.581) . we may think that f/e„ {T) —r Uq (T) weakly-star in X°°, for some Uq G X°° 
such that Uq satisfies 

\\Uo (T)llxoo < Q(l + II /ill I LOO (f^) + ||/i2|lLoo(r)). (3.59) 

This together with (I3.56P gives U (T) G X°° such that U (T) satisfies the analogue of (I3.59p . 

Finally, it remains to prove that U (T) is also bounded in by essentially arguing as in 
the proof of Proposition 13. 7[ Note once again that in contrast to the limit case e„ = 0, the 
parabolic system (I2.17P is non-degenerate if e„ > 0, and we have enough regularity of Ue„ 
to justify the multiplication by the test functions dtUe„ in the weak formulation (12. lip (cf. 
Definition EH). Analogous to fl335D - fl338D . we get 

dtS,^ (t) + 2 \\dtu,„ (t)||J,(^) + 2 \\dtv,„ ml2^r,ds/b) = 0' (3-60) 

for the energy S^^ defined in (I3.36p . Therefore, Gronwall's inequality applied to (I3.60p yields 
as in dSSHD-dSSi, 

(T) + (\\dtu,^ (r)lli.(^) + \\dtv,^ {r)\\l.^r4S/b)) dr (3.61) 

< Q{1 + ||/j-i||ioo(j^) + ||^2|lLoo(r)), 

which in light of (I3.37P gives that 11^^ (T) is uniformly (in T, e„) bounded in V-". Passing now 
to the limit as e„ — > 0, we obtain in a standard way that U (T) satisfies (I3.6ip as well. The 
desired inequality (I3.55P follows immediately from (I3.59P and (I3.6ip . The last assertion in 
the theorem is a standard corollary of the energy identity fl2.38p . So, the proof is complete. 
■ 

The (forward) uniqueness theorem holds for bounded solutions of (ll.ip -( ITT2|) . (I2.12p . 

Theorem 3.17. Let the assumptions of Theorem \3.1(A he satisfied, and let f,g & (M, M). 
Consider two functions Ui (t) , U2 [t) G C ([0, T] ; X^) which solve / fO]) -/ f7r^) . i\2.12\) in the 
sense of Definition \2Jl In addition, let Ui,U2 G L~ ([0, T] ; X~). Then, f/i (0) = f/2 (0) 
implies that Ui (t) = U2 (t) , for all t G [0, T] . 

Proof. The proof is an immediate consequence of inequality ^2.^5^ (cf. Proposition \2.^) . 
Indeed, for E [t) := \\Ui (t) — U2 (t) H^a , we have E E (0, T), and 

\dtE (t) < Q (||f/i (t)||^^ , \\U2 mW^) E (t) < cE (t) , (3.62) 
which yields the desired claim on account of the application of Gronwall's inequality. ■ 
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Remark 3.18. We note that in contrast to the non-degenerate case of a (s) = v, we do not 
know whether the backwards uniqueness theorem holds for hounded solutions of the parabolic 
system hl.l\) - [T7^) forp ^ 2. Namely, if the equality Ui (T) = U2 (T) holds for some T >0, 
then we have Ui (t) = U2 (t) , for t < T, as well. Indeed, in the former case the parabolic 
system U.l\) - l[T^) is just a reaction- diffusion equation with dynamic boundary conditions, so 
we can establish additional regularity of the weak solutions in L°° Y^''^)nW^'°° (M+; V^'^), 
follo wing f2^ . Theorem 2.3] (see also below). Thus, exploiting a well-known theorem (see. 



e.g., 14 B . Theorem 11.10]; cf. also JiB, Chapter HI]), we can easily establish the backwards 



uniqueness result in this case. 

Finally, we observe that since a bounded weak solution U (t) G X°° fl V^) is 

unique, we may define a global attractor Agi for the parabolic system fll.ip -f lTT^ . f l2.12p by 
the standard expression 

Agi := HoAr, (3.63) 
and define a classical semigroup on this attractor via 

St : Agi ^ Agi, StU (0) = f/ (t) , (3.64) 

Here, U (t) is the unique (bounded) weak solution of fll.ip - fll.2p . fl2.12p . such that U (t) 
satisfies the energy identity fl2.38p . We also note that estimate ( 13.550 gives a uniform estimate 
of the L°°-norm of the trajectories belonging to the attractor Agi. Therefore, the growth rate 
of the nonlinearities /, g with respect to U becomes nonessential for further investigations 



of global attractors and we can study them exactly as in, e.g., [42|, |45 



3.3. Strong trajectory attractors: the semilinear case 

Let (T2 (t) , /C_|_) be the trajectory dynamical system associated with the reaction-diffusion 
equation fll.4p . subject to the dynamic boundary condition 01.50 (see Section 1). In this 
section, we shall establish additional regularity estimates for the weak solutions of (11.41) - 
(II. 5p . and obtain as a by-product, that the weak trajectory attractor Atr, constructed in 
Theorem I3.14[ is in fact a strong trajectory attractor. In order to do so we will verify, for 
every U G Atr, that the attraction property and the compactness holds not only in the 
weak topology of Q^'''°'^, but also in the strong topology of G^'"'^. The definition of a strong 
trajectory attractor is obtained by replacing the weak attraction condition (iii) in Definition 
13.131 by the condition of strong attraction in the topology of G^'°^. 

We have the following proposition, whose proof goes essentially as in jis]. Theorem 2.3]. 

Proposition 3.19. Let the assumptions of either Theorem \3.S\ or Proposition \3.4\ be satis- 
fied, and f,g & (R, R) . Every bounded complete weak solution U & K, of problem jl.j^ - 
[TR) belongs to (R; V^'^) n W^'"^ (R; V^'^) , and the following inequality holds for t e R, 

11^7 (t) 11^2,2 + \\dtU{t)\\l,., + j^^ \\dlU{r)\\l,dr (3.65) 

< Q{1 + ||/;.i||ioo(j^) + ||^2||Loo(r)), 
where Q is a monotone nondecreasing function, independent of t, U and the data. 
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Proof. Indeed, having established the L°°-estimate fl3.55p . fl3.65p can be easily derived 
using a standard technique for parabolic equations with dynamic boundary conditions (see 
23, Theorem 2.3] for further details; cf. also 10, 25, 2^|). ■ 
Consequently, we have shown the following. 

Theorem 3.20. Under the assumptions of Proposition \3A^ the weak attractor Atr, con- 
structed in Theorem 3. 14, is a strong trajectory attractor for the trajectory dynamical system 
(T2 (t) , /C+). Moreover, Atr is compact in Cioc V^"'''^) , for any s G (0, 1], and the set 

dtAtr ■■= {dtU : U e Atr} 
is compact in Cioc V-*^^''^) , for any I G (0, 1/2). 

Proof. It is sufficient to show that the set T (t^) i3, where > 1 is sufficiently large and 
B is the absorbing set f l3.5ip for the semigroup {T (t)} is compact in the strong topology of 
the space 0^''°^^. We observe that for every t > t# and any U G /C+, we have 

T{t)U e L°° (M+; V^'^ n X~) n W^'°° (M+; V^'^) n (M+; X^) . (3.66) 

Thus, it is enough to prove that the kernel /C+, for the attractor Atr defined by f l3.50p . is 
compact in the strong local topology of the space O^. We first note that, due to estimates 
( 13351) . (Km and ( Km . we have 

/C+ is bounded in Cb (M+; V^''^'^) , for any s G (0, 1], 
and moreover, it is compact in the local topology 

}C+cQoc{^+;^'-'''). (3.67) 
This follows from the following embedding 

{U:UeLZ V^'^) , dtU G L~ (M+; V^'^) } c Q^c V^-^'^) . 

which is compact. Thus, in view of fl3.67p and the boundedness of /C+ in L°° (M+; X°°), we 
immediately see that /C+ is also compact in the (strong) local topology of 0']_, i.e., we have 

/C+ C e^;,,,, := (M+; X^) n LL (M+; V^'^) n {LZ, (M+; {Q)) x (M+; L^' (r))) . 
The second statement follows analogously using the compactness of the following embedding 

{dtU : U G (M+; V^'^) , d^U G LL (M+; X^) } c V^-''^) . 

The proof is finished. ■ 

Corollary 3.21. Under the validity of assumptions of Theorem \3.2(j\ the reaction- diffusion 
equation (T^, with the dynamic boundary condition U.5\) possesses a global attractor Agi{= 
/C (0)), defined by ^M^-^EM' ^^^^ A« bounded m V^'^ n X°°. 



37 



Remark 3.22. It is worth mentioning that one can also establish more regularity of the weak 
solution u E (fi) as much as it is allowed by the regularity of VL, f, g and the external 
sources hi, /i2. Taking advantage of the regularity result in Corollary VJ. 21\ we can prove that 
the global attractor Agi is finite dimensional, by establishing the existence of a more refined 
object called exponential attractor £gi. However, since the associated solution semigroup St 
happens to be (uniformly quasi-) differentiable with respect to the initial data, on the attractor 
Agi, we can instead employ a volume contraction argument (see, e.g., l23ll ). 



We can now extend the results in |23| for the case of nonhnear boundary conditions, 
without requiring that the restrictive condition (H2) holds. More precisely, we can establish 
the following upper-bound on the dimension of the global attractor Agi. 

Theorem 3.23. Provided that f,g E C^(M, M) satisfy the assumptions of Corollary \3.21\ 
for as long as ^3.2\) holds, the fractal dimension of Agi = /C (0) admits the estimate 

dimp Agi < Ci (1 + ly-^^-^^) , forN>2 (3.68) 

and 

dimp Agi < Ci (l + z/-^/2) ^ forN = 1. (3.69) 
where Ci depends only on Q, T and the sources hi,h2. 

Proof. First, it is easy to establish that the flow St : Agi — )■ Agi generated by the reaction- 
diffusion equation (11. 4p and dynamic boundary condition (II. 5p is uniformly differentiable on 
Agi, with differential 

L (t, U (t)) : $ = (^^^^ G ^ 1/ = Q G X^, (3.70) 

where V is the unique solution to 

dtv = uAv - f {u{t))v, (^dt'^ + iybdnv + g'{u{t)^^)v^^^ = 0, (3.71) 
V{0) = $. 

Indeed, the uniform differentiability result follows from the assumptions on /, g and is a 
consequence of the boundedness of Agi into V^'^ fl X°° (see (iij). In order to deduce (I3.68P - 
(I3.69p . it is sufficient (see, e.g., [ij. Chapter III, Definition 4.1]) to estimate the j-trace of 
the operator 

iyA-f'{u{t)) 



L(t,f/(t)) - , , „ , , , ^, 

' ' -budn -g {v{t)) 
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We have 



Trace (L {t, U (t)) Q„,) = ^1 it, U {t)) y^,-, 



X2 



2,r 



m 



J = l 



^ (/' (m [t]) ^j, ifj) - ^ {v (t)) 



2,r 



where the set of real- valued functions ipj G fl V^'^ is an orthonormal basis in Q^X^. 
By Theorem 13.161 it follows that every bounded complete trajectory U{t), t G M, for the 
dynamical system (5^, K, (0)) is uniformly bounded in V^'^ fl X°°, namely, it holds: 



sup \\U 



< C 



(3.72) 



where the positive constant C is independent of U (t) , u, Q, T, but depends on the L°°-norms 
of the external forces hi, /i2 and the constants in (13.71) (this follows from the usual description 
of the global attractor, see, e.g., [oj, and the main estimates). Thus, exploiting (12.361) once 
more there holds 



max sup 



L°°(n) 



g' iv{t)) 



L°=(r) 



for some > which depends on C and the growth rates ri, r2 > 1; we find 



(3.73) 



Trace {L {t,U) Qm) < -uj^ 11^^^3111 + 



m. 



From |23|, Proposition 5.5], we obtain 



Trace (L [t, U) Qm) < -vciCw T) m^+^ + [ciuCw T) + CJm 
=: p{m) , 

for some absolute positive constant ci independent of the parameters of the problem {Cw 
is given explicitly in 23|, Theorem 5.4]). The function p{y) is concave. The root of the 
equation p{y) = is 

/ r \^-^ 

Thus, we can apply a well-known result [14, Chapter VIII, Theorem 3.1] to deduce that 
dim p Agi < y*, from which (I3.68P follows. The case = 1 is similar. ■ 
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Remark 3.24. In the competing scenario l[2. we can also explicitly estimate from above 
the dimension of the global attractor Agi = K, (0) for the dynamical system {St,Agi) defined 
by /i3.6!^) - f3.64\ )- However, in this case the upper bound does not seem to be as sharp as in 



i\3. 68\) - l[3. 69\) . Indeed, by Theorem \3.16[ (see also Theorem \3.2\) . it follows that every bounded 



complete trajectory U{t), t G M, for the dynamical system {St, )C{0)) satisfies ( [5*. 72\ ) with 



constant C = ^ v"^ as v ^ 0^, for some (Xi < 0; by assumption (H3), this implies that 



the constant in \3. 13^ behaves as = C^^iy) ~ v'^'^ , for some 02 < —1 depending on 
(Ti and ri,r2. Hence, in this case there seems to be a discrepancy between the upper bound 
( [5'. 74\ ) and the lower bound in U.13\) . 



3.4- A blowup result 

As pointed out at the beginning of this article, nonhnear dissipative boundary conditions 
cannot prevent blowup of some solutions of fll.4l) - fll.5p when the non-dissipative interior term 
/ is superlinear, i.e., when / satisfies (11.91) for some ri > 2. We will follow some arguments 



similar to ones presented in [43| for nonlinear Robin boundary conditions, by constructing 
some subsolutions to (ll.4p - (ll.5l) which become unbounded in finite time at some points of 
the boundary F. We begin with the following notion. 

Definition 3.25. A function v : Q x [0,T) ^ W is a subsolution of llil.4\ )- fl~5\) if it satisfies 

dtv -uAv + f (v) - hi (x) < 0, mnx (0, T), 
dtv + ubdr,v + giv)- h2 (x) < 0, onT x (0, T) ^ ' 

and 

V (0) < Mo in n. (3.76) 



Analogously, the function v is called a supersolution if the inequalities in ^3.75^ - [3. 76\ ) are 
reversed. 



From [4, Section 7], we have the following 

Proposition 3.26. Let u be a solution of ( [j.^[ j- pT3]) . and let v and v be a subsolution and 
supersolution, respectively, of ( [j.^| )-( (i.5]) . in the sense of Definition \3. 25i Then, 

V (x, t) <u (x, t) <v (x, t) , 

for all X E fl and for as long as they exist. In particular, if / (0) < and g (0) < 0, and if 
Uq > 0, then the solution of ( [j./^| )- pT3]) satisfies u {x,t) > 0, for all x E Vl, for as long as it 
exists. 



We aim to construct subsolutions by comparing solutions of (ll.4l) - (ll.5p with classical 
solutions fulfilling certain Dirichlet conditions on the time lateral boundary F x (0,T). For 
that purpose, the following result is very useful (see, 
follows, it suffices to consider the case hi = 0, h2 = 0. 



cf. also 0, Eij). In what 
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Proposition 3.27. Assume that there are a C^-concave decreasing function h{s) and a 
number Sq > such that 

limsup/i (s) < —Xi^D, 

s— J-oo 

where Xi^d > is the first eigenvalue of —uA with Dirichlet boundary conditions. Moreover, 
suppose that h{s) < for all s > sq such that 

< oo. (3.77) 



so 



\h(s 



Then, there are positive (smooth and locally well-defined) solutions v of the reaction- diffusion 
equation 

dtv - uAv + f{v) = 0, inVLx (0, T), (3.78) 
subject to the boundary and initial conditions 

V = Q, onV X (0,T) , ,^ 

V (0) := fo, inVL, 

that blowup in finite time. 

Theorem 3.28. Let the assumptions of Proposition \J7El\ be satisfied, and assume in addition 
that f (s) < h (s) < for all s > sq such that ( [J. y?| j holds. Then, for any nonlinear function 
g, there exist solutions of ( [i.^p - pT3]) that blowup in finite time. Moreover, there exists a 
positive function Wq (x) such that all solutions with initial data Uq greater or equal than 
Wq + vq, blowup in finite time. 

Proof. Let ipi be the principal eigenfunction associated with Xi d > such that ||v'i|lii(n) = 
1. It is well-known that > in by the maximum principle. Thus, by choosing A : = 
max {so, Sq, O} such that h (s) < —Xi d, for all s > Sg, we can define w (x) := 5(fi (x) > 0, 
for some 6 > 0, such that 

i^bdnW + ^ (A) < on r. (3.80) 

Note that it is always possible to fix 5 > since, by the maximum principle once again, we 
have dn'^i < on F. Let us now define 

U (x, t) = Wq (x) + V (x, t) , 

where Wq (x) := w (x) + A and f is a solution of fl3.78p - fl3.79p . Arguing as in the proof of 



43l . Theorem 4.3], we can easily establish that 

dtu - uAu + f (u) <0, infl X (0, T) . 
On the other hand, on F x (0,T), we have 

dfU + ubdnU + g (u) = vbd^w + vbd^v + g {A) 
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but since d-^v < on F, by the choice of A and w (cf. fl3.80p ). it follows 

dtu + vhdryU + c/ (m) < 0, on r X (0, T) . 

Consequently, we deduce that m is a subsolution of f ll.4p -( IT3|) . Therefore, on account of 
Propositions 13.26 1 and I3.27t all solutions u = u {x,t) of fll.4l) - fll.5p that satisfy Uq > Wq + Vq 
in Q, blowup in finite time. The proof is complete. ■ 

Corollary 3.29. Assume that f satisfies 

lis) 
s (In s) 

for some I > 1. Then the conclusion of Theorem \3. 2S\ applies. 



limSUp — j < 0, 



4. Appendix 

In this section, we state two basic results which are essential to the analysis of problem 
fll.ip - fll.2l) . The first lemma is just a variation of a result in 28|, Section 5] using well-known 
facts about nonlinear forms and maximal monotone operators in Sobolev spaces. For the 
convenience of the reader, we give below a simple proof of that result. 

Lemma 4.1. Let a (s) satisfy all the assumptions in Theorem \2.2\ and let U be fixed. 
Then the functional V Bp{U,V), defined by l{2.13\) . belongs to (V^)*. Moreover, Bp is 
strictly monotone, hemicontinuous and coercive. 



Proof. We will make use of 28|, Definitions 2.1 and 2.2]. In our case, for p > pq, 



V:=YP CX^= (X^)* cV* = {YPy . 

VP be fixe( 

fl2.3p . we obtain 



Let U = (^^) G VP be fixed. It is clear that Bp{U, ■) is hnear. Let V = {Jj G V^. Exploiting 



This implies that Bp{U, ■) G (V^)* , for every U G V^. 
Next, let U,V e Y^. Then, recaUing ([23D, we have 



BpiU,V)\<\\ur^'j\V\W.. (4.1) 



Bp{U,U-V)-Bp{V,U-V) (4.2) 

= / {a{\Vuf)Vu-a{\Vvf)Vv) ■V{u-v)dx+ I {\u\''-'^u ~ \v\''-\) {u - v)dx 
Jn Jn 

> / {\u\ + \v\f~^\u-v\'^dx > 0, 
Jn 



which shows that Bp is monotone. This estimate also shows that 



BpiU,U-V)-BpiV,U-V)>0, 
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for all U,V G with U ^ V. Thus, Bp is strictly monotone. The continuity of the norm 
function implies that Bp is hemicontinuous. Finally, it is easy to deduce that 

Ita = ,4.3) 

||{/||yp— >+oo 



which shows that Bp is coercive. The proof is finished. ■ 

The following lemma is useful in controlling surface integrals by means of volume integrals. 
The proof follows from an application of the divergence theorem to the function w = {\u\^ 
n, for a smooth vector field E (Q, M^) such that ^ • n = 1. 

Lemma 4.2. Let p > 1, s > 1 and u G W^'^ (Q) . Then for every e > 0, there exists a 
positive constant = C (e, , independent of u, such that 

where 7 = max (s — 1) / (p — 1)) . 
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